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PnLL-aURFACB  IKTBRFBROMBTRIC  SCJ^HMER 
FOR  THB  MBABURBMBHT  OP  QRAZIHQ  INCIDBKGB  MIRRORB 


1.0  IMTRODUGTIOH 

The  objective  of  the  Phase-II  research  was  to  design,  construct  and 
test  a  prototype  interferometer  system  (FSIS)  for  the  purpose  of 
testing  cylindrical-like  mirrors  such  as  those  used  in  grazing 
incidence  X-ray  imaging  applications. 

Existing  instruments  are  based  on  a  shearing  interferometry,  dual 
laser  beam  scanning  approach,  measuring  the  slope  or  curvature 
functions  along  a  single  long  trace  in  the  long  direction  of  the 
cylindrical  mirror.  While  these  Instruments  are  characterized  by  a 
high  measuring  accuracy,  they  are  very  slow  and  environmentally 
sensitive.  An  even  slower  and  more  complex  two-dimensional  scan 
must  be  implemented  if  the  full  surface  measurement  is  required. 

Our  approach  is  also  based  on  interferometric  slope  measurement  in 
the  direction  of  the  long  side  of  the  mirror,  using  a  grating 
shearing  interferometer  head.  A  full-width  beam  covering  an 
extended  subaperture  of  the  mirror  is  scanned  in  discrete  steps  in 
the  direction  of  the  long  side  of  the  mirror,  so  that  a  single  scan 
allows  coverage  of  the  whole  mirror  surface.  By  least  square 
fitting  the  series  of  subaperture  measurements,  the  system  yields 
parallel  measurement  of  a  multiplioity  of  long  traces  in  the 
direction  of  scan. 

In  the  following  final  report  we  describe  the  work  on  the  prototype 
FSIS  system  including  a  description  of 

The  optical  interferometer  subsystem 

The  digital  and  electronics  hardware 

The  data  reduction  software. 

We  then  describe  the  testing  program  and 
results  for  a  cylindrical  test  surface. 
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2.0  PRINCIPLES  AND  ANALYSIS 


In  classical  Twyman-Green  or  Fizeau  type  interferometry,  the 
wavefront  under  test  is  compared  with  an  ideal  wavefront  reflected 
from  a  reference  plane  or  spherical  surface.  Systems  based  on  these 
approaches  are  capable  of  yielding  very  high  measuring  accuracy, 
but  are  extremely  sensitive  to  environmental  disturbances  such  as 
vibration  and  turbulence.  Moreover,  they  are  characterized  by  a 
fixed  sensitivity  and  thus  cannot  deal  with  strongly  aberrated 
surfaces  that  result  in  high  fringe  densities. 

Our  FSIS  system  makes  use  of  lateral  shearing  interferometry  using 
a  coarse  Ronchi  grating  as  the  shearing  device.  In  this  approach, 
there  is  no  need  for  a  reference  surface  since  the  wavefront  under 
test  is  made  to  interfere  with  a  laterally  shifted  replica  of 
itself.  The  wavefront  under  test  is  passed  through  the  grating  so 
that  interference  is  obtained  between  the  overlapping  +1  and  -1 
diffraction  orders  produced  by  the  grating.  The  result  is  an 
interferometer  that  is  greatly  insensitive  to  vibration  and  whose 
sensitivity  can  be  varied  by  changing  the  pitch  of  the  grating. 
However,  the  measured  guantity  is  now  the  slope  of  the  wavefront 
rather  the  wavefront  itself  as  in  classical  interferometry  and  an 
additional  numerical  integration  step  is  reguired  to  obtain  the 
desired  wavefront  function. 

Choosing  a  coordinate  system  where  x  and  y  are  respectively  in  the 
directions  of  the  short  and  long  side  of  the  cylindrical-type 
mirror,  and  orienting  the  grating  with  its  pitch  in  the  direction 
of  the  long  side  of  the  cylindrical-type  mirror,  the  inter ferogram 
intensity  is  given  by  the  following  expression: 

I(x,y)  =  A  +  B  Cos  [Phi(x,  y  s/2)  -  Phi(x,  y  -  s/2)]  , 

where  s  is  the  shear  value  and  Phi  is  the  phase  of  the  wavefront 
under  test. 

The  measured  quantity  is  thus  equal  to 

Phi(x,  y  +  s/2)  -  Phi(x,  y  -  s/2)  =  s  (dPhi/dy)  , 

i.e.  the  slope  of  the  wavefront  in  the  y  direction  as  desired.  The 
wavefront  Phi(x,y)  is  then  obtained  by  numerical  integration  of  the 
slope  data 

Phi(x,y)  =  1/s  [Phi(x,  y  +  s/2)  -  Phi(x,  y  -  s/2)]  dy  . 

A  measurement  of  a  long  cylindrical-type  mirror  is  obtained  by 
scanning  the  mirror  with  an  full-width  cylindrical  wavefront  over 
a  finite  number  of  overlapping  subapertures.  The  full-surface 
wavefront  is  then  obtained  by  numerical  synthesis  by  least  square 
fitting  of  the  series  of  subaperture  measurements.  The  result  of 
the  measurement  is  a  collection  of  wavefront  profiles  or  traces  in 
the  direction  of  the  long  side  of  the  mirror. 


-2- 


3.0  THE  OPTICAL  IMTBRFEROMBTgR  BUBBY8TBM 

3.1  BvateM  Operation 


A  schematic  diagram  of  the  optical  interferometer  system  is  shown 
in  Figure  l  and  includes  the  following: 

-  A  coherent  laser  source 

A  spatial  filter  to  'clean  up*  and  expand  the  beam,  Including 
a  microscope  objective  mounted  on  an  X-Y  translation  stage  and 
focused  onto  a  rotating  diffuser  plate 

A  lens  LI  to  collimate  the  laser  beam  diverging  from  the 
spatial  filter 

A  lens  L2  to  focus  the  beam  onto  the  grating  shearing  device 

A  coarse  Ronchi  grating  having  a  couple  of  lines  per 
millimeter,  mounted  on  an  X-Y-Z  translation  stage  and  a 
rotation  stage  for  focusing  and  alignment  purposes 

A  high  quality,  large  collimating  lens  L  to  produces  plane 
wave  illumination 

A  high  quality  cylindrical  lens  to  produce  the  desired 
cylindrical  wavefront  used  to  illuminate,  in  a  null 
configuration,  the  cylindrical  type  surface  under  test  at 
normal  incidence  as  shown  in  top  and  side  views 

A  beam  splitter  plate  that  recombines  the  reflected  beam  onto 
the  CCD  array  detector  through  some  Imaging  optics 

A  High  resolution  CCD  camera  to  capture  the  interferograms. 

The  cylindrical  mirror  under  test  is  mounted  onto  a  long  travel 
linear  motorized  stage.  The  stage  is  translated  in  discrete  steps 
so  that  the  mirror  is  fully  covered  by  a  number  of  overlapping 
subapertures  as  shown  in  Figure  2. 

The  Y-axis,  of  the  grating  translation  stage  is  equipped  with  a 
micro-stepper  stage  that  permits  to  translate  the  grating  in  its 
plane  to  vary  the  phase  of  the  interferogram  by  small  fractional 
increments . 

Both  motorized  stages,  as  well  as  the  CCD  camera  are  under  computer 
control  so  that  the  operations  of  scanning  and  data  acquisition  are 
Implemented  in  a  fully  automated  fashion. 
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3.2  The  Optical  Prototype 


Based  on  the  design  of  Figure  1,  an  integrated  prototype  system 
was  designed  and  constructed. 

A  diagram  of  the  prototype  system  is  shown  in  Figure  3,  including: 

The  laser  source:  2mW  helium-neon  laser  operating  at  632. 8nm 

The  spatial  filter  with  high  power  objective  and  rotating 
diffuser 

Beam  shaping  optics:  collimating  lens  LI  and  focusing  lens  L2 
assembly 

Beam  splitter  plate 
Grating  shearing  device 

Reference  beam  forming  optics  for  null  testing:  Collimating 
lens  L  and  cylindrical  lens  CL  assembly 

Imaging  lens  L3 

-  CCD  camera  detector. 


Figures  4  to  7  show  photographs  of  the  prototype  system: 

Figure  4  shows  a  photograph  of  the  complete  optical  head  and 
computer  control 

Figure  5  shows  the  microstepper  stage  (ORIEL)  which  controls 
the  grating  position 

Figure  6  shows  the  spatial  filter  assembly  with  the  rotating 
diffuser  unit 

Figure  7  shows  a  close-up  photograph  of  the  long  travel 
translation  stage  (NEAT)  with  the  cylindrical  mirror  under 
test  mounted  via  a  tilt/rotation  stage  for  alignment  purposes. 
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Figure  3:  Schematic  Diagram  of  Interferometric  Head 
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Figure  Photograph  of  Complete  FSIS  Prototype  System 


Figure  5;  Photograph  of  FSIS  System  Showing  the  Oriel  Stepper 
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FiiiURE  7:  Photograph  of  the  FSIS  System  Showing  the  Test  Mirror 
Assembly 
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4.0 


DIGITAL  AMD  ELECTRONICS  HARDWMIB 


The  digital  and  electronics  hardware  for  the  FSIS  includes 

The  control  electronics  hardware:  motor  controllers  to  drive 
the  test  mirror  scanning  stage  and  the  grating  microstepper 
stage 

The  data  acquisition  hardware:  CCD  camera  and  frame  buffer  to 
acquire,  digitize  and  store  the  inter ferograms 

The  data  processing  and  display  hardware:  microcomputer  system 
to  process  the  inter ferograms  to  yield  surface  data  and 
display  the  results. 


Diagram  1  shows  a  schematic  of  the  complete  digital  and  electronics 
hardware  used,  including: 

Micr  jcomputer 

Hard  and  floppy  disk  drives 

-  VGA  monitor 
Keyboard  and  mouse 

-  Laser  printer 
Video  printer 

Motor  controllers  (NEAT  and  ORIEL) 

Frame  grabber  board 
CCD  camera . 
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Diagram  1:  Block  Diagram  of  the  Digital/Electronics  System 
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4.1  Control  Blectronlca  Hardware 


The  control  electronics  hardware  includes  the  motor  controllers  to 
drive 


The  long  travel  scanning  stage  (NEAT)  that  holds  and  translate 
the  mirror  under  test 

The  short  travel  microstepper  (ORIEL)  that  shifts  the  grating 
shearing  device  by  fine  incremental  steps. 


4.1.1  Scanning  Stage  Control  fNEATl 

The  NEAT  translation  stage  is  a  long  travel  (24  Inches)  motorized 
stage  under  computer  control  as  shown  in  Diagram  2  below. 

The  detailed  specifications  of  the  controller  and  the  wiring 
connections  with  the  computer  are  described '  in  the  following 
section  entitled  ‘NEAT  Translation  Stage  Control  System*. 


4.1.2  Microstepper  Stage  Control  System  (ORIELS 

The  Oriel  Microstepper  stage  is  used  to  shift  the  grating  shearing 
device  within  the  interferometer  head.  It  is  controlled  through  the 
computer  serial  port  RS232  COMl,  as  shown  in  Diagram  3. 

The  detailed  specifications  of  the  controller  and  the  wiring 
connections  with  the  computer  are  described  in  the  following 
section  entitled  'ORIEL  Microstepper  Control  System* . 
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BLOCK  DIAGRAM  OP  THE  NEAT  BCAMHINQ  STAGE  CONTROL 


STEPPIMQ  MOTOR  CONTROLLER  SYSTEM 


Diagram  2 .  Block  Diagram  of  the  NEAT  Scanning  Stage  Control  System 


[EAT  TRANBIATIOM  BTAQB  COMTROL  SYSTEM 

NEAT-310M  Controller  Electronlc/Mechanloal  Specification 

Drive  Type: 

54  volt  bipolar  chopper 
.7  to  3.5  Amps/Phase,  half-coil 
1.5  to  7  Amps/Phase,  full-coil. 

Memory  buffer: 

32  K  bytes  non-volatile  program  memory 
Velocity  Range: 

40  to  327,640  steps  /  sec. 

Position  Range: 

-8,388,608  to  -*-8,388,608  steps 
Computer  interface: 

RS  232  and  8-bit  parallel  communication  ports 
Power  requirements: 

115/230  V  A.C.  50/60  Hz,  1.5  Amps 


MEAT-310M  Controller  CQM«iiii<qatlon  Wiring  Format 

Computer:  RS-232  9  pin  (Female)  adapter. 

Pin  #2:  (Red) 

Pin  #3:  (Green) 

Pin  #5:  (Black) 

Pin  #6,7,8:  Short  circuit. 

Controller:  25  pin  (Male)  adapter. 

Pin  #2:  (Green) 

Pin  #3:  (Red 
Pin  #1,7:  (Black) 
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3. 


HEAT-3 lOM  Software  CQMiHwnil 


NEAT-3 lOM  controller  uses  RS-232  serial  communications  port  to 
handle  control  commands  from  the  PC  computer.  The  communication 
port  setting  are:  8  data  bit,  1  stop  bit,  no  parity  bit,  ASCII  data 
format. 

The  NEAT-3 lOM  operates  with  a  command  set  of  over  60  high-level 
ASCII  commands.  The  following  examples  are  common  command  sets  we 
are  using  in  the  control  software  package: 

MR:  Move  relative. 

Format:  MRddddd  ddddd  range  from  -8,388,608  to  +  8,388,608 
MA:  Move  absolute. 

Format:  MAddddd  ddddd  range  from  -8,388,608  to  +  8,388,608 

MH:  Move  physical  origin  position. 

MC:  Move  continuous. 

Format  MC+/-. 

SP:  Set  absolute  register  counter  value. 

VI:  Set  initial  velocity.* 

Format:  VIddddd  ddddd  range  from  40  to  327,680 
default  is  4000  steps  /  sec. 

-  VF:  Set  final  velocity. 

Format:  VFdddddddd  dddddddd  range  from  40  to  327,680 

AC:  Set  acceleration  and  deceleraticii. 

Format:  AC (+/-) dddddddd  ,  *+*  for  acceleration,  for 

deceleration,  without  sign  for  both. 

ST:  Stop  moving  immediately. 

-  ME:  Mnemonic  Expansion. 

Format:  MFE  for  enable,  MED  for  disable. 

MF:  Move  finish  signal  feed  back. 

Format:  MFE  for  enable,  MFD  for  disable. 


*  Note:  The  primary  resonance  is  at  400-1200  steps  /  sec.  the 
maximum  allowable  start/stop  frequency  is  load  dependent,  typically 
below  6000  /  sec  is  safe. 
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BLOCK  DIAGRAM  OF  THE  ORIEL  MICROSTEPPER  COMTROL 


8TBPPIHQ  MOTOR  CONTROLLER  SYSTEM 


Diagram  3.  Block  Diagram  of  ORIEL  Microstepper  Control  System 
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ORIBL  MICROSTBPPER  STAGE  CONTROL  SYSTEM 


1.  Controller  Bleotronlo /Mechanical  Specif Icatlon: 

Number  of  controlled  stepper  motors: 

Up  to  two,  individually. 

Data  I/O  : 

300,600,1200,2400,4800,9600  baud,  switch  selectable. 
Compatible  stepper  motor: 

2  or  4  phase  unipolar  motor,  24V,  0.5  A  max,  per  phase. 
Winding  holding  current: 

Reduced  to  0.1  A  per  phase . 

Limit  switch  connections: 

Low  logic  level  =  limit  reached. 

Maximum  speed: 

Half  step  mode:  1000  steps  /  sec. 

Full  step  mode:  500  steps  /  sec. 

Computer  interface: 

RS-  232  serial  communication  port 

Power  requirements: 

110/220-240  VAC  50/60  Hz,  switch  selectable. 

2 .  Stepper  Mlfce  Motor 

Model:  ORIEL  Stepper  Mike.  Model  18500. 

Specification: 

Step  Size: 

Half  step  mode:  1.3,  0.7  urn. 

Full  step  mode:  2  urn. 

Maximum  step  rate: 

Half  step  mode:  1000  steps  /  sec. 

Full  step  mode:  500  steps  /  sec. 
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Maximum  spindle  speed: 

1  mm  /sec. 

Maximum  axial  load: 

15.5  lbs. 

Uni-directional  repeatability: 

<  2  urn. 

Backlash: 

<  3  urn. 

Travel  Range: 

0.5”,  (13  mm) . 


3.  Oriel  ateppar  Mike  Controller  cation  Port  Wiring  Format 

Computer:  RS232  9  pin  (F)  adapter. 

Pin  #2:  (Red) 

Pin  #3:  (Green) 

Pin  #5:  (Black) 

Pin  #6,7,8:  short  circuit. 

-  Controller:  25  pin  (M)  adapter. 

Pin  #2:  (Red) 

Pin  #3:  (Green) 

Pin  #1,7:  (Black) 

Controller  to  Motor: 

10  lines  Rainbow  cable. 


4.  Oriel  Stepper  Mike  Controller  Boftware  cnmand  Bets 

Oriel  Stepper  Mike  Controller  uses  PS-232  serial  communications 
port  to  handle  control  commands  from  the  PC  computer.  The 
communication  port  setting  are:  7  data  bit,  1  stop  bit,  no  parity 
bit,  ASCII  data  format. 

The  controller  operates  with  command  sets  simular  to  NEAT-3 lOM.  The 
following  list  are  some  of  the  common  command  sets  we  will  use  in 
the  control  oftware  package. 
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A,B;  Direct  the  flow  of  cominand,  data,  status  Inquiries  to 
motor  control  register. 

C:  Clear  absolute  register  to  zero,  (set  logic  origin) . 

D:  Disable  motor  driver. 

E:  Enable  the  motor  driver. 

F:  Change  to  full  step  mode. 

G:  Go  absolute. 

Format  G  +/-  dddd,  is  necessary  to  stop  data  field,  dddd 

is  the  number  of  steps  to  go. 

H:  Change  to  half  step  mode. 

I:  Inquire  the  status  of  motor,  return  'd',  d  =  0  -  7 

bit  0,  direction:  CW  =  1  ,  CCW  =  0. 

bit  1,  step  size:  half  =  1  ,  full  =  0. 

bit  2,  E/D  status:  Enable  =  1,  disable  =  0. 

Q:  Inquire  the  status  of  translator,  return  'd',  d  =  0  -  7 
bit  0,  motor:  On  =  l  ,  Off  =  0. 

bit  1,  CW  limit:  yes  =1  ,  no  =  0. 

bit  2,  CCW  limit:  yes  =1  ,  no  =  0. 

R:  Set  up  the  desired  step  rate  for  each  motor. 

Format:  Rn,  n  *  0  -  7. 

S:  To  start  the  moving. 

T:  Travel  relative  steps. 

Format :  Tdddddddd , 

Z:  Move  to  the  logic  origin. (the  position  of  absolute 
register^  0) 

§:  Unconditional  stop  both  motor. 

>:  Move  single  step  in  CW  direction. 

<:  Move  single  step  in  CCW  direction. 

=:  Start  both  motors  simultaneously. 

*  After  command  set,  eg  AT1000,BT1000. 

+:  Set  direction  to  CW. 

-:  Set  direction  to  CCW. 
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4.2  Data  Acquisition  Hardware 


The  data  acquisition  hardware  includes: 

A  high-resolution  imaging  CCD  camera  for  interferograro 
detection 

A  frame  grabber  and  buffer  for  A/D  conversion  and  storage 

A  TV  monitor  for  real-time  interferogram  viewing  and  a  video 
printer  for  obtaining  interferogram  hard  copying. 

A  block  diagram  of  the  complete  data  acquisition  system  is  shown 
in  Diagram  4. 

4.2.1  The  CCD  Camera 

The  CCD  camera  used  for  interferogram  data  capture  is  an  area 
camera  with  RS-170  video  signal  output  made  by  COHU,  INC. 

The  specifications  of  the  camera  are  as  follows: 

Imager:  Single  CCD  using  the  frame  transfer  method 

-  Image  Area:  8.8mm  x  6.6mm  or  2/3  format 

Active  Pixels:  754  x  484 

Resolution:  565  horizontal  x  350  vertical  TV  lines 

Sensitivity:  0.07  lux  useable  picture, 

0.16  lux  at  full  video 

Signal/Noise:  56  dB 

Contrast  Variation:  Less  than  5%  overall  @  25  degree  C 
Gamma:  Jumper  selectable  at  1  or  o.45 

AGC:  Jumper  selectable,  on/off 

signal:  EIA  RS-170  standard. 
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4.2.2  The  Frame  Grabber /Buffer  Board 


The  frame  grabber/buffer  board  used  as  an  interface  for  the  CCD 
camera  is  the  Silicon  Video  Mux  made  by  EPIX,  INC. .  The  Silicon 
Video  Mux  is  a  single  board  frame  grabber  that  allows  the  PC 
computer  to  digitize,  process,  display,  transmit  and  archive  video 
information. 

The  board  can  digitize  one  or  several  frames  from  a  video  camera 
such  as  a  CCD  camera  and  allow  the  PC  to  process  the  image  data, 
and  display  the  image  data  on  a  black/white  or  RGB  monitor.  Video 
data  is  digitized  and  displayed  at  8  bits  per  pixel,  and  the  board 
can  contain  up  to  4  MB  of  image  memory  for  multiple  image  storage. 

The  specifications  of  the  frame  buffer  board  are  as  follows: 

Variable  sampling 

1  MB  or  4MB  image  memory 

Video  memory  address  registers 

752  pixels  by  480  pixels  maximum  resolution 

External  pixel  clock,  horizontal  and  vertical  drive 

-  Multiple  image  storage 

Video  line  counters  for  selection  of  raster  lines  to  be 
digitized/displayed. 


4.2.3  Output  Devices 

The  Interferograms  are  viewed  in  real  time  on  a  TV  monitor 
connected  directly  to  the  CCD  camera  or  through  the  frame  buffer. 

A  color  or  black/white  video  printer,  connected  to  the  monitor,  is 
used  to  obtain  hard  copies  of  the  interferograms. 

A  video  recorder,  also  be  connected  to  the  monitor,  can  be  used  to 
record  sequences  of  inter ferogram  data. 
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CONNECTION  DIAGRAM  OF  VIDEO  INPUT/OUTPUT  SYSTEM 
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COLOR 

VIDEO  PRINTER 

TV  MONITOR 

VIDEO  RECORDER 

Diagram  4. 


Block  diagram  of  the  video  input/output  system 


4.3  Data  Process ing  Hardware 


An  IBM  compatible  microcomputer  system  is  used  for  interferogram 
data  processing. 

The  system  Includes: 

-  Computer  with 

-  80486/40  MHs  Microprocessor 

64K  Cache  Memory 
8  MB  RAM 

Hard  Disk  Drive:  200  MB 
Dual  Floppy  Drives: 

5.25*'  Drive  /  1.2  MB 
3.5"  Drive  /  1.44  MB 

-  Ports : 

1  Parallel  Port 

2  Serial  Ports 
Display  Card: 

VGA 

Memory:  1  MB 
Monitor: 

VGA  Color 

0.28mm  Non-Interlaced  with  1024x768  Resolution 
101  Enhanced  Keyboard 

-  Microsoft  Mouse 
MS-DOS  5 

Laser  Printer:  HP  LaserJet  III. 
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5.0 


SOFTWARE 


5.1  SOFTWARE  OVERVIEW 

The  Multitrace  system  software  package  Includes  the  following 
files: 

0.  Installation  program:  It  Installs  the  package  software  Into 
the  IBM  PC  compatible  computer. 

PROGRAM :  INSTALL . BAT 

1.  Batch  file  for  compiling  main  program:  It  allows  the  user  to 

modify  the  main  program  package  and  recompile  Into  the 

executable  file. 

PROGRAM :  COMPILEl . BAT 

2.  Main  program  of  the  package:  This  program  controls  the 

Multitrace  Profiler  system,  acgulres  the  test  object  data, 
process  and  display  the  results. 

PROGRAM :  PACKAGE . BAS 

PROGRAM :  PACKAGE . EXE 

3.  Batch  file  for  compiling  demo  program:  It  allows  the  user  to 

modify  the  demo  program  package  and  recompile  Into  the 

executable  file. 

PROGRAM :  COMPI LE2 . BAT 

4.  Demo  program  of  the  package:  This  demo  program  controls  the 
Multitrace  Profiler  system,  process  the  pre-stored  the  test 
object  data  and  display  the  results. 

PROGRAM :  DEMO . BAS 

PROGRAM :  DEMO . EXE 

5.  Data  set  for  demo  purpose. 

PROGRAM :  DEMOI . DAT 

PROGRAM :  DEMOR . DAT 

6.  Image  acgulsltlon  subroutine,  (called  by  main  program) 

PROGRAM:  ACQ.C 

PROGRAM:  ACQ.EXE 

7.  Compiler  support  libraries  and  subroutines. 


PROGRAM :  BRUN30 . LIB 

PROGRAM :  USERLIB . EXE 
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5.2  User  Manual :  Operation  of  tha  Multi-trace  Svatem 


The  cylindrical  mirror  testing  package  user  manual  includes 
the  hardware  motion  controllers  specification  sheets  with 
sample  programs  to  test  the  motion  control  and  system 
operation  menu  driven  software  package  with  the  programs. 

The  motion  control  specifications  list  the  model  and  maker  of 
the  controllers  and  motors,  their  electrical  and  mechanical 
specifications,  wire  connections  of  the  communication  ports 
and  control  software  command  lists. 

The  system  operation  package  is  menu  driven  type  of  functions. 
When  type  "MAINF”  in  DOS  system,  the  display  will  show  a 
message  and  the  MAIN  MENU  on  the  screen.  By  pressing  the 
Function  keys  for  the  desired  operation,  the  user  can 
manipulate  the  interferometer  system  and  implement  some 
specific  control  and  data  processing  for  the  test  mirrors. 
There  are  two  operations  in  the  main  menu  and  are  Implemented 
by  pressing  function  keys  FI,  F2,  pressing  the  ''ESC  key  will 
exit  the  program  and  back  to  DOS  system. 

The  MAIN  MENU  is  shown  in  the  following: 


-  FSI8  SYSTEM 

FI  :  TEST  MIRROR  SCAMNXNa 
F2  I  FULL  SURFACE  MEASUREMENT 
ESC:  EXIT 


The  functions  of  each  operation  are  explained  in  the 
following: 


FI.  TEST  MIRROR  SCAMMINQ 

The  TEST  MIRROR  SCANNING  is  to  operate  the  interferometer 
system  in  the  scanning  mode,  which  allows  the  user  to  align 
the  cylindrical  test  mirror  with  the  Interferometer  and  scan 
the  full  surface  of  the  test  mirror  to  view  the  complete 
inter ferogr am. 

When  Function  key  FI  is  pressed,  the  linear  translation  stage 
(NEAT)  will  be  activated  to  move  the  test  mirror  along 
X~direction  (cylindrical  mirror  axis)  by  five  subaperture 
steps  back  and  forth  at  a  certain  speed  to  allow  the  user  to 
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view  and  simulate  the  image  acquisition  of  the  full 
interferogram  of  the  test  mirror.  This  processing  can  be 
repeated  as  many  times  as  wish  by  pressing  the  FI  as  needed. 


F2.  FULL  SURFACE  MEASUREMENT 

The  FULL  SURFACE  MEASUREMENT  is  to  operate  the  interferometer 
system  to  scan  the  cylindrical  test  mirror  in  a  multi-aperture 
mode.  For  the  first  test  mirror,  it  is  assigned  five 
subaperture  to  cover  the  full  mirror. 

It  takes  the  first  subaperture  inter ferograms  and  stores  in 
the  frame  buffer  then  move  the  test  mirror  to  the  2nd 
position,  then  takes  the  2nd  subaperture  inter ferograms  and 
stores  in  the  frame  buffer.  The  same  procedure  repeats  until 
all  five  subaperture  interferograms  are  acquired  and  stored  in 
the  frame  buffer. 

After  all  the  subaperture  interferograms  being  acquired.  The 
computer  starts  to  process  the  phase  maps  of  each  subaperture 
and  synthesize  these  subaperture  phase  into  full  surface  phase 
map.  Then  the  screen  display  the  3-D  plot  of  the  synthesis 
slope  phase  map  and  allow  the  user  to  key  in  the  section 
number  on  the  surface  and  plot  the  section  trace  profile. 

To  escape  from  the  Interactive  sectional  profile  screen,  just 
type  '•999”  follow  the  "ENTER"  to  move  to  next  display  "THIN 
LINE  CONTOUR".  Press  "ENTER"  to  leave  this  page. 

The  same  display  and  Interactive  procedure  for  synthesis 
wavefront  of  the  test  mirror. 


ESC:  MAIM  MENU 

When  "ESC"  key  is  pressed,  the  program  will  return  back  to 
MAIN  MENU  page. 
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5.3 


SOFTWARE  PACKAGE  flOURCB  CODE 


ITSKI  'h  **  h  It  It  It  ■k'k'k  hit  *11  it  It  it  It  hit  It  hit 

rem  *  INSTALL.BAT  * 

IT  sin  'klfkltltltltltltliltlfkltlfkltlflilfitlfk 

@echo  off 
echo  Caution  ! ! ! 

echo  To  prevent  from  version  conflict,  this  installation  will 
echo  delete  all  pre  exist  files  echo  in  the  destination  directory, 
echo  Ctrl-Break  to  stop  installation  ;  Or 
pause 

if  '%i'==:''  goto  nodir 
c:  >  NUL 
md  %1  >  NUL 
cd  %1 

echo  Delete  pre-exist  files 
del  <  yes  >  NUL 

echo  installing. . . 
copy  b:\*.*  >  NUL 

goto  done 
: nodir 
c:  >  NUL 

md  \multra  >  NUL 
cd  \multra 

echo  Delete  pre-exist  files 
del  *.*  <  yes  >  NUL 
echo  installing... 
copy  b:\*.*  >  NUL 

:done 

echo  Installation  is  done, 
echo  on 


rem  **************..  •%*★***** 
rem  *  C0MPILE1.BAT  * 
rem  aaaaaaaaaaaaaaaaaaaaaaa 

qb  package  /I  /d  /q  ; 
link  package  ; 


rem  *  COMPILE2.BAT  * 

rem  *********************** 

qb  demo  /I  /d  /q  ; 
link  demo  ; 


•  ***************************** 

’  *  PACKAGE. BAS  * 

*  ***************************** 
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'This  Is  the  complete  program  package  for  the  multitrace  system 


COMMON  SHARED  selkey  'declare  command  variable 
const 
thold%=200 ,  cutr%=10 ,  cutc%=20,  ton%=l,  tof  f  center%=120,  slze%-240 

, lnltord%=0 

const  scalepl=. 5729583, smth%=l 
const  shlft%  =  3 

const  stepr%='3,  xscale=4.2,yscale=2.5,yshlft=3.8,stpno%  =  3 

const  xdc  =  20,ydc  =  170,111%  =  l,rrr%  =  120 

common  shared  top% , bot% , lb% , rb% 

common  shared  max,mln,buf% 

common  shared  row%,set% 

common  shared  coel,coe2,coe3 

common  shared  maxp,mlnp, maxi, mini, rmsp,rmsl,pvp,pvl,rms,pv 
common  shared  levellt, zscale,deglnc,mark%, tltle$ 
rem  $dynamlc 

dim  shared 

Obj (240) ,c%(240) ,r%(240) ,ZZZ(120,5) ,tz(360) ,ttz(360) , shlftpts% (5) 

dim  shared 

order%(240) , slope (120, 120) ,zz (120, 120) ,xb%(2,120) ,bb%(2,5) , jff%(5 

) , jee%(5) ,fllenm$(5) 

dim  shared  colls%(8) 

dim  shared  m(480) ,mp(480) 

dim  shared 

z(150) ,h(9) ,lsh(9) ,xh(5) ,yh(5) ,lm(9) ,x(120) ,y(120) , jm(5) ,castab(3 
.3,3) 

rem  $statlc 

d  a  t  a 


0,0, 8, 0,0, 2, 5, 1,7, 6, 9, 2, 3, 4, 5, 6, 0,3, 4, 7, 1,0, 1,8, 4, 3, 0,9, 0,1, 2, 3, 9 
,6, 7, 4, 5, 2, 0,5, 8, 0,0, 6, 7, 8, 9, 1,3, 5, 0,3, 8, 0,0, 3, 8, 0,3, 6, 9, 0,1, 1,0 


screen  9 


call  prlntmenu 
CALL  Initial 
DO  UNTIL  selkey  =27 
program 

CALL  kbhlt 
operator 

SELECT  CASE  selkey 
CASE  59 
to  each 

CALL  send("MR+18750") 
scanning 

call  stopped 

CALL  send("MR+18750*') 

call  stopped 

CALL  8end(''MR+18750") 

call  stopped 

CALL  send("MR+18750") 

call  stopped 

CALL  send  (•'MR-75000'’) 

call  stopped 


'ESC  key  to  exit 
'wait  for  key  hit  from 

'have  the  stage  moved 
'aperture  for  testing 
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call  printmenu 
CASE  60 
CLS 

LOCATE  12,10 

print  "  ACQUISITION  SUBAPERTURE  #1 

N 

shellCacq  5”)  'call  acq.exe  to  get  the 

image 

call  waiting 
CALL  sendCMR+lSTSO") 
call  stopped 
LOCATE  12  10 

print  "  "  ACQUISITION  SUBAPERTURE  #2 

N 

ShellCacq  10") 
call  waiting 
CALL  send("MR+18750") 
call  stopped 
LOCATE  12,10 

print  "  ACQUISITION  SUBAPERTURE  #3 

M 

ShellCacq  15") 
call  waiting 
CALL  send("MR+18750") 
call  stopped 
LOCATE  12,10 

print  "  ACQUISITION  SUBAPERTURE  #4 

N 

shell("acq  20") 
call  waiting 
CALL  send("MR+18750") 
call  stopped 
LOCATE  12,10 

print  "  ACQUISITION  SUBAPERTURE  #5 

If 

shell ("acq  2.  ' 
call  waiting 
CALL  send ("MR-, 5000") 
call  stopped 
call  main 
call  printmenu 
CASE  ELSE 
END  SELECT 
LOOP 

cls 

END 

•  ** 

•***  Dummy  waiting  subroutine. 

•  ** 

sub  waiting  static 

for  il  «  -32767  to  32767 
for  jj  =  0  to  1 
next  jj 
next  li 
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end  sub 


I  ** 

•***  main  process,  contains  computation  ,  and  display. 

•  ** 

sub  main  static 
call  init 
call  onoff(toff%) 
for  sett  ■  1  to  5 
buft  ■  sett  *  5 
max  *  -32767 
min  «  32767 
locate  12,20 

print  "  COMPUTING  SUBAPERTURE'';sett;" 

n 

call  compute 
next  sett 
call  onoff(tont) 
call  dodata 
call  display 
end  sub 


•  ** 

'***  wait  for  asterisk  to  make  sure  command  is  accept  by  controller 

•  ** 

SUB  ASTERISK  STATIC 

DO  UNTIL  INPUT$(1,  06)  *=  "*•' 

LOOP 

END  SUB 


•  ** 

'***  initial  process 
•  ** 

SUB  initial  STATIC 

OPEN  "COM2;9600,N,8,1,RS,CS,DS,CD'»  FOR  RANDOM  AS  #6  'OPEN  SERIAL 
PORT 

CALL  send ("MED")  'disable  mnemonics 

CALL  send("MFE")  'enable  move  finish 

CALL  send("VI600")  'set  initial  velocity 

CALL  send("VF6000")  'set  final  velocity 

CALL  send("AC1001")  'set  acceleration  and 

deceleration 

END  SUB 

•  ** 

•***  accept  key  hit  from  keyboard 
•  ** 

SUB  kbhit  STATIC 
key$  -  "" 

DO 
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WHILE  key$  =  "” 
key$  -  INKEY$ 

WEND 

IF  LEN(key$)  >  1  THEN 

skey$  »  MID$(key$,  2,  1) 

ELSE 

skey$  »  MID$(key$,  1,  1) 

END  IF 

selkey  »  ASC(skey$) 

LOOP  UNTIL  key$  <> 

END  SUB 

•  ** 

•***  subroutine  print  the  operation  menu  onto  screen 

•  ** 

SUB  prlntmenu  STATIC 
CLS 

LOCATE  3,1 

print  "  CYLINDRICAL  MIRROR  MULTITRACE  MEASUREMENT 

N 

print  "" 
print  "*• 
print 
print 

print  •*  FI  :  TEST  MIRROR  SCANNING” 

print 

print  "  F2  :  ACQUISITION  AND  MEASUREMENT” 

print  ”” 
print  "” 

print  ”  ESC:  EXIT” 

END  SUB 

•  ** 

•***  subroutine  for  sending  command  to  controller 

•  ** 

SUB  send  (a$)  STATIC 

PRINT  #6,  a$  'send  command  to 

controller 

CALL  ASTERISK  'wait  until  asterisk 

is  return 
END  SUB 

•  ** 

'***  wait  until  process  finish 
•  ** 

SUB  STOPPED  STATIC 

DO  UNTIL  INPUT$(1,  #6)  «  ”F” 

LOOP 

for  li  »  -32767  to  32767 
for  jj  “  0  to  1 
next  jj 
next  il 

END  SUB 


-34- 


•  ** 

•  *** 
•  ** 

sub 


initialize  parameter, 

init  static 
filenm$(l)  =*  "rO.dat" 
=  "rl.dat" 
»  "r2.dat •• 
*  "r3.dat" 
»  ••r4.dat" 
=  10 
=  3 

*  9 

=  1 
»=  4 

*  12 
=  14 
*  15 
im(l)=l 
jm(l)=0 


end 


filenm$(2) 

£llenm$(3) 

£llenm$(4) 
f llenm$ (5) 
colls% (1) 
colls% (2) 
colls% (3) 
colis%(4) 
colls%(5) 
colisi (6) 
colis%(7) 
colls%(8) 
im(0)=0  ;  im(l)=l  :  im(2)=l 
jm(0)=0  :  jm(l)=0  :  jm(2)=l 
for  k=0  to  3 
for  j=0  to  3 

for  i*0  to  3 
read  castab(k, j , i) 
next  i 
next  j 
next  k 
sub 


im ( 3 ) =0 
jm(3)»l 


•  ** 

'***  initialize  parameter  of  slope  display, 
sub  initp  static 

title$="CYLINDRICAL  MIRROR  :  MULTITRACE  MEASUREMENT  /  SLOPE 
FUNCTION" 

levell%  *  12 
zscale  =  -.3 
deginc  »  36 
mark%  -  60 
for  1%  =  1  to  480 
m(i%)  »  32767 
mp(i%)  =  32767 
next  i% 
max  ■■  maxp 
min  «  minp 
rms  «  rmsp 
pv  «  maxp  -  minp 
end  sub 


•  ** 

<***  initialize  parameter  of  wavefront  display. 

•  ** 

sub  initi  static 

title$-"CYLINDRICAL  MIRROR  :  MULTITRACE  MEASUREMENT  /  WAVEFRONT 
FUNCTION" 
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levell%  =  180 
zscale  *  -.005 
deglnc  =  720 
inark%  »  1800 
for  1%  =  1  to  480 
in(i%)  »  32767 
inp(i%)  =  32767 
next  1% 
max  »  maxi 
min  -  mini 
rms  =  rmsi 
pv  *  maxi  -  mini 
end  sub 


'***  compute  slope  of  each  aperture. 

•  ** 

sub  compute  static 

open  fllenm$(set%)  for  output  as  #7 
call  baseorder 
if  top%  <  cutr%  *  2  then 
top%  =  cutr% 
else 

top%  =  top%  /  2  +  1 
end  if 

if  bot%  >  (120  -  cutr%)  *  2  then 
bot%  =  120  -  cutr% 

*bot%  =  bot%  /  2  -  1 
end  if 

for  l%=top%  to  bot% 
row%  =  i%  *  2  +  1 
call  getrowphase 
if  lb%  <  cutc%  *  2  then 
xb%(l,l%)  =  cutc% 

6lS6 

xb%(l,i%)  *  lb%  /  2  +  1 
end  if 

if  rb%  >  (120  -  cutc%)  *  2  then 
xb%(2,i%)  =  120  -  cutc% 

6lS6 

xb%(2,i%)  =  rb%  /  2  -  1 
end  if 

for  j%  **  xb%(l,l%)  to  xb%(2,i%) 

if  j%  *  2  >  lb%  and  j%  *  2  <  rb%  then 
if  obj(j%  *  2)  <>  -32767  then 
zz(j%,i%)  *  obj(j%*2) 
else 

zz(j%,i%)  »  -32767 
end  if 
else 

zz(j%,i%)  =  -32767 
end  if 
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next  j% 
next  i% 
call  smooth 
call  rotate 

print  #7, top%,bot%,max,mln 
for  i%  «=  top%  to  bot% 
print  #7,l%,xb%(l,l%) ,xb%(2,l%) 
for  j%  =  xb%(l,l%)  to  xb%(2,l%} 
print  #7,slope( j%, 1%) 
next  j% 
next  1% 

close 
end  sub 

•  ** 

•***  trace  the  order  of  phase. 

•  ** 

sub  traceord  static 

order % ( center % ) =lnltord% 
cc%=lnitord% 

for  j%*center%  to  slze%-l 
If  obj ( j%+l)-obj ( j%)<-180.0  then  cc%=cc%+l 
If  obj (j%+l) -obj (j%) >180.0  then  cc%=cc%-l 
order % ( j  %+l ) =cc% 
next  j% 
cc%*lnltord% 

for  j%*center%  to  1  step  -1 
If  j%<>center%  then 

If  obj (j%+l) -obj (j%) <-180.0  then  cc%=cc%-l 
If  obj(j%+l)-obj(j%)>180.0  then  cc%=cc%+l 
end  If 

order%( j%)*cc% 
next  j% 
end  sub 

•  ** 

•***  convert  the  phase. 

'  ** 

sub  phaseconverslon(t%(l) )  static 
for  j%=l  to  240 
If  t%(j%)  <>  -32767  then 

ol^j  ( j*)*t%  ( j%)  *scalepl  /  2 
else 

obj(j%)  »  -32767 
end  If 
next  j% 
end  sub 

•  ** 

I***  get  the  base  order. 

*** 

sub  baseorder  static 

call  ptr86(se9a%,  offa%,  varptr(c% (1) ) ) 

call 

getpc ( top% , bot% , thold% , slze% , slze% , buf % , center% , sega% , of f a% ) 
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call  phaseconverslon(c%() ) 
call  traceord 
end  sub 

•  ** 

'***  get  phase. 

•  ** 

sub  getrowphase  static 

call  ptr86(sega%,  offa%,  varptr (r% (1) ) ) 

call 
getpir ( lb% , rb% , thold% , order! ( row!) , size! , size! , buf ! , row! , sega! , of 
fa!) 

call  phaseconverslon(r!() ) 
end  sub 

•  ** 

•***  data  smoothing. 

•  ** 

sub  smooth  static 

for  i!  =  1  to  120 
for  j!  =  1  to  120 

Slope(j!,l!)  =  -32767 
next  j! 
next  i! 

for  i!  «=  top!  to  bot! 
for  j!=xb!(l,i!)  to  xb!(2,i!) 
if  zz(j!,l!)  <>  -32767  then 
sum  =  0 
cnt  =  0 

for  ii!  =  i!  -  smth!  to  i!  +  smth! 
if  ii!  >  top!  and  ii!  <=  bot!  then 
for  jj!  =  j!  -  smth!  to  j!  +  smth! 
if  jj!  >  xb!(l,ii!)  and  jj!  <  xb!(2,li!)  then 
if  zz(jj!,ii!)  <>  -32767  then 
sum  =  sum  +  zz(jj!,ii!) 
cnt  =  cnt  +  1 
end  if 
end  if 
next  jj! 
end  if 
next  ii! 

slope(120-i!, j!)  =  sum  /  cnt 
if  max  <  slope (120-i!,j!)  then 
max  *  slope (120-i!,j!) 
else 
end  if 

if  min  >  slope(120-l!, j!)  then 
min  «  slope(120-i!, j!) 
else 
end  if 
end  if 
next  j! 
next  1! 
end  sub 
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•  ** 

•***  data  rotating. 

•  ** 


sub  rotate  static 
ntop%  =  32767 
nbot%  =  -32767 
for  i%  =  top%  to  bot% 

if  xb%(l,i%)  <  ntop%  then  ntop%  =  xb%(l,i%) 
if  xb%(2,i%)  >  nbot%  then  nbot%  =  xb%(2,l%) 
next  i% 

for  i%  =  1  to  120 
xb%(l,i%)  =  121 
xb%(2,l%)  =  0 
next  i% 

for  i%  =»  ntop%  to  nbot% 
for  j%  =  1  to  120 

if  slope(j%,i%)  <>  -32767  then 
xb%(l,i%)  =  j% 
j%  =  121 
end  if 
next  j% 

for  j%  =  120  to  1  step  -1 

if  slope(j%,i%)  <>  -32767  then 
xb%(2,i%)  =  j% 
j%  =  0 
end  if 
next  j% 
next  1% 

for  1%  =  ntop%  to  nbot% 
if  xb%(l,l%)  <  xb%(2,i%)  then 
top%  =  i% 
i%  =  121 
end  if 
next  i% 

for  i%  =  nbot%  to  ntop%  step  -1 
if  xb%(l,l%)  <  xb%(2,i%)  then 
bot%  *  i% 
i%  =  0 
end  if 
next  i% 
end  sub 


•  ** 

'***  match  data  to  a  complete  window. 
•  ** 

sub  dodata  static 
for  i%  =  1  to  5 

shlftpts%(i%)  =  shift!  *  (i%  -  1) 
next  1% 

max  *  -32767  :  min  =  32767 
open  "rO.dat"  for  input  as  il 
open  "rl.dat"  for  input  as  #2 
open  ''r2.dat”  for  input  as  #3 
open  ”r3.dat"  for  input  as  #4 
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open  '’r4.dat''  for  input  as  #5 

input  #l,topl%,botl%, maxi, mini 
input  #2,top2%,bot2%,max2,min2 
input  #3 , top3%,bot3%,max3 ,min3 
input  #4,top4%,bot4%,max4,min4 
input  #5,top5%,bot5%,max5,min5 
top%  =  0 


if  topl% 

> 

top% 

then 

top% 

= 

topl% 

if  top2% 

> 

top% 

then 

top% 

ss 

top2% 

if  top3% 

> 

top% 

then 

top% 

= 

top3% 

if  top4% 

> 

top% 

then 

top% 

= 

top4% 

if  top5%  > 
bot%  =121 

top% 

then 

top% 

— 

top5% 

if  botl% 

< 

bot% 

then 

bot% 

= 

botl% 

if  bot2% 

< 

bot% 

then 

bot% 

= 

bot2% 

if  bot3% 

< 

bot% 

then 

bot% 

= 

bot3% 

if  bot4% 

< 

bot% 

then 

bot% 

= 

bot4% 

if  bot5% 

< 

bot% 

then 

bot% 

= 

bot5% 

for  i%  =  topl%  to  top%  -  1 
input  #l,dmy%,bl%,b2% 
for  j%  =  bl%  to  b2% 
input  #l,z% 
next  j% 
next  i% 

for  i%  =  top2%  to  top%  -  1 
input  #2,dmy%,bl%,b2% 
for  j%  *  bl%  to  b2% 
input  #2,z% 
next  j% 
next  i% 

for  i%  =  top3%  to  top%  -  1 
input  #3,dmy%,bl%,b2% 
for  j%  =  bl%  to  b2% 
input  #3,z% 
next  j% 
next  i% 

for  i%  =  top4%  to  top%  -  1 
input  #4,dmy%,bl%,b2% 
for  j%  =  bl%  to  b2% 
input  #4,z% 
next  j% 
next  i% 

for  i%  =  top5%  to  top%  -  1 
input  #5,dmy%,bl%,b2% 
for  j%  =  bl%  to  b2% 
input  #5,z% 
next  j% 
next  i% 

for  i%  ®  top%  to  bot% 

input  #l,dmy%,bb%(l,l) ,bb%(2,l} 
nrb%  =  120  -  cut 
for  j  '  «  bb%(l,l)  to  bb%(2,l) 
if  j*  >  nrb%  then 
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input  #l,dmy% 
else 

Input  #l,zzz(j%,l) 
end  If 
next  j% 

if  bb%(2,l)  >  nrb%  then  bb%(2,l)  =  nrb% 

input  #2,diny%,bb%(l,2)  ,bb%(2,2) 
for  j%  =  bb%(l,2)  +  shiftpts%(2)  to  bb%(2,2) 
if  j%  <  1  or  j%  >  120  then 
input  #2,diny% 
else 

input  #2,zzz(j%,2) 
end  if 
next  j% 

if  bb%(l,2)  <  1  +  cut  then 

bb%(l,2)  =  cut  +  1  +  shiftpts%(2) 
else 

bb%(l,2)  =  bb%(l,2)  +  shiftpts%(2) 
end  if 

if  bb%(2,2)  >  120  -  cut  then 

bb%(2,2)  =  120  -  cut  +  shiftpts%(2) 
else 

bb%(2,2)  =  bb%(2,2)  +  shlftpts%(2) 
end  if 

if  bb%(2,2)  >  120  then  bb%(2,2)  =  120 


input  #3,diny%,bb%(l, 3)  ,bb%(2,3) 
for  j%  =  bb%(l,3)  +  shiftpts%(3)  to  bb%(2,3) 
ifj%<lorj%>  120  then 
input  #3,droy% 
else 

input  #3,zzz(j%,3) 
end  if 
next  j% 

if  bb%(l,3)  <  1  +  cut  then 

bb%(l,3)  =  cut  +  1  +  shiftpts%(3) 
else 

bb%(l,3)  =  bb%(l,3)  +  shlftpts%(3) 
end  if 

if  bb%(2,3)  >  120  -  cut  then 

bb%(2,3)  =  120  -  cut  +  shiftpts%(3) 
else 

bb%(2,3)  =  bb%(2,3)  +  shiftpts%(3) 
end  if 

if  bb%(2,3)  >  120  then  bb%(2,3)  =  120 


input  #4,diny%,bb%(l,4)  ,bb%(2,4) 
for  j%  *  bb%(l,4)  +  shiftpts%(4)  to  bb%(2,4) 
if  j%  <  1  or  j%  >  120  then 
input  #4,diny% 
else 


+  shiftpts%(2) 


+  shlftpts%(3) 


+  shiftpts%(4) 
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Input  #4,zzz(j%,4) 
end  If 
next  j% 

if  bb%(l,4)  <  1  +  cut  then 

bb%(l,4)  *  cut  +  1  +  shlftpts%(4) 
else 


bb%(l,4) 
end  if 
if  bb%(2,4) 
bb%(2,4) 
else 
bb%(2,4) 
end  if 
if  bb%(2,4) 


=  bb%(l,4)  +  shiftpts%(4) 

>  120  -  cut  then 

=  120  -  cut  +  shiftpts%(4) 

=  bb%(2,4)  +  shlftpts%(4) 

>  120  then  bb%(2,4)  =  120 


input  #5,dmy%,bb%(l,5) ,bb%(2,5) 
for  j%  =  bb%(l,5)  +  shiftpts%(5)  to  bb%(2,5) 
if  j%  <  1  or  j%  >  120  then 
input  #5,diny% 
else 

input  #5,zzz(j%,5) 
end  if 
next  j% 

if  bb%(l,5)  <  1  +  cut  then 

bb%(l,5)  =  cut  +  1  +  shiftpts%(5) 
else 


bb%(l,5)  =  bb%(l,5)  +  shiftpts%(5) 
end  if 

if  bb%(2,5)  >  120  -  cut  then 

bb%(2,5)  *  120  -  cut  +  shiftpts%(5) 
else 


bb%(2,5)  *  bb%(2,5)  +  shlftpts%(5) 
end  if 

if  bb%(2,5)  >  120  then  bb%(2,5)  =  120 
for  ii%  =  1  to  5 
for  j%  =  1  to  bb%(l,ii%)  -  i 
zzz(j%,ll%}  =  -32767 
next  j% 

for  j%  =  bb%(2,ii%)  +  1  to  120 
zzz(j%,ii%)  =  -32767 
next  j% 
next  ii% 


for  j%  =  1  to  360 
tz(j%)  =  -32767 
next  j% 

for  ii%  »  1  to  4 
ni%  =  il%  +  1 
sum  “  0 
cnt  *  0 

jff%(ii%)  *  bb%(l,ni%)  +  60 

jee%(ii%)  =  bb%(2,ii%) 

for  j%  -  jff%(il%)  to  jee%(ii%) 


+  shiftpts%(5) 
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nj%  =  j%  -  60 

if  2Z2(j%,ii%)  <>  -32767  and  222(nj%,ni%)  <>  -32767  then 
sum  =  sum  +  (22z(j%,ii%)  -  222 (nj%,nl%) ) 
cnt  =  cnt  +  1 
end  if 
next  j% 

if  cnt  >  0  then 
dev  -  sum  /  cnt 
else 
dev  =  0 
end  if 

for  j%  =*  bb%(l,ni%)  to  bb%(2,ni%) 

if  222(j%,ni%)  <>  -32767  then  222(j%,ni%)  =  222(3%, nit)  + 

dev 

next  j% 
pit  =  iit  -  1 
jdct  =  pit  *  60 
if  iit  =  1  then 

for  j%  =  bb%(l,ii%)  to  jff%(iit) 
t2(jt)  =  Z22(jt,iit) 

•print  iit, jt,t2 ( jt) 
next  jt 

else 

for  jt  =  jeet(pi%)  -  60  to  jff%(ii%) 
t2(jt+jdct)  *  222(jt,ii%) 

•print  iit, j%+jdc%,t2(jt+jdc%) 
next  jt 
end  if 

for  jt  =  jff%(ii%)  to  jee%(iit) 
njt  =  jt  -  60 

if  (222(j%,ii%)  <>  -32767  and  222(nj%,ni%)  <>  -32767)  then 
t2(j%+jdct)  =  (222(jt,iit)  +  222(nj%,ni%) )  /  2 
else 

if  222 (jt, iit)  <>  -32767  then 
t2(jt+jdct)  =  222(jt,iit) 
else 

if  222(njt,nit)  <>  -32767  then 
t2(jt+jdc%)  =  222(njt,nit) 
else 

t2(jt+jdct)  =  -32767 
end  if 
end  if 
end  if 
next  jt 

if  iit  =  4  then 

for  jt  *  jee%(ii%)  -  60  to  bbt(2,ni%) 
t2 ( jt+jdc%+60)  =  222(jt,nit) 

•print  iit, j%+jdct+60,t2( jt+jdct+60) 
next  jt 
end  if 
next  iit 
2one%  »  3 
smptst  ®  5 
ptpt  *  2 
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ptn%  *  3 

for  ii%  *  1  to  4 
sum  «  0 
cnt  =  0 

ff%  *  jff%(ii%)  +  60  *  (11%  -  1)  -  zone% 

ee%  =  jff%(ll%)  +  60  *  (11%  -  1)  +  zone% 

for  j%  =  ff%  -  ptp%  to  ff%  +  ptn%  -  1 
If  tz(j%)  <>  -32767  then 
sum  =  sum  +  tz(j%) 

cnt  =  cnt  +  1 

end  If 
next  j% 

for  j%  =  ff%  to  ee% 

If  tz(j%)  <>  -32767  then  ttz(j%)  =  sum  /  cnt 
If  tz(j%  -  ptp%)  <>  -32767  then 
sum  =  sum  -  tz(j%  -  ptp%) 
cnt  =  cnt  -  1 
end  If 

If  tz(j%  +  ptn%)  <>  -32767  then 
sum  =  sum  +  tz(j%  +  ptn%) 
cnt  *  cnt  +  1 
end  If 
next  j% 

for  j%  =  ff%  to  ee% 

If  tz(j%)  <>  -32767  then  tz(j%)  *=  ttz(j%) 
next  j% 
sum  «  0 
cnt  =  0 

ff%  -  jee%(ll%)  +  60  *  (11%  -  1)  -  zone% 

ee%  =  jee%(ll%)  +  60  *  (11%  -  1)  +  zone% 

for  j%  =  ff%  -  ptp%  to  ff%  +  ptn%  -  1 
If  tz(j%)  <>  -32767  then 
sum  =  sum  +  tz(j%) 

cnt  =  cnt  +  1 

end  If 
next  j% 

for  j%  =  ff%  to  ee% 

If  tz(j%)  <>  -32767  then  ttz(j%)  =  sum  /  cnt 
If  tz(j%  -  ptp%)  <>  -32767  then 
sum  =  sum  -  tz(j%  -  ptp%) 
cnt  =  cnt  -  1 
end  If 

If  tz(j%  +  ptn%)  <>  -32767  then 
sum  =  sum  +  tz(j%  +  ptn%) 
cnt  =  cnt  +  1 
end  If 
next  j% 

for  j%  =  ff%  to  ee% 

If  tz(j%)  <>  -32767  then  tz(j%)  =  ttz(j%) 
next  j% 
next  11% 

for  j%  »  1  to  120 

zz(j%,l%)  =  -32767 
next  j% 
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for  j%  =  4  to  117 

zz(j%,l%)  =  tz(j%  *  3) 

if  zz(j%,i%)  <>  -32767  then 

if  max  <  zz(j%,i%)  then  max  =  zz(j%,i%) 

if  min  >  zz(j%,i%)  then  min  =  zz(j%,i%) 

end  if 

next  j% 

xb%(l,i%)  =  121 
xb%(2,i%)  =  0 
for  j%  =  1  to  120 

if  zz(j%,i%)  <>  -32767  then 

xb%(l,i%)  =  j% 

j%  =  121 

end  if 

next  j% 

for  j%  =  120  to  1  step  -1 

if  zz(j%,l%)  <>  -32767  then 
xb%(2,l%)  =  j% 
j%  =  0 
end  if 
next  j% 
next  i% 
call  tiltterm 
call  removetilt 
call  integrate 
end  sub 


'***  integrate  the  wavefront. 

t  ** 

sub  integrate  static 
maxi  =  -32767 
mini  *=  32767 
sumi  =  0  :  cnti  =  0 
for  i%  *  1  to  120 
for  j%  «  1  to  120 

zz(j%,i%)  =  -32767 
next  j% 
next  i% 

for  i%  *  top%  to  bot% 
zz(60,i%)  =  0 
for  j%  =®  61  to  xb%(2,l%) 

if  slope(j%,i%)  <>  -32767  then 

zz(j%,i%)  =  zz(j%-l,i%)+slope(j%,i%)  /  1.5 
6lS6 

zz(j%,i%)  =  zz(j%-l,i%) 
end  if 

if  maxi  <  zz(j%,i%)  then  maxi  =  zz(j%,i%) 
if  mini  >  zz(j%,i%)  then  mini  =  zz(j%,i%) 
sumi  =  sumi  +  zz(j%,i%)''2 
cnti  “  cnti  +  1 
next  j% 

for  j%  =  59  to  xb%(l,i%)  step  -  1 
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if  slope (j%,i%)  <>  -32767  then 

Z2(j%,i%)  =  2z(j%+l,i%)-slope(j%+l,i%)  /  1.5 
else 

zz(j%,l%)  =  zz(j%+l,i%) 
end  If 

if  maxi  <  zz(jt,i%)  then  maxi  =  zz(j%,i%) 
if  mini  >  zz(j%,i%)  then  mini  =  zz(j%,i%) 
sumi  =  sumi  +  zz(j%,i%)^2 
cnti  =  cnti  +  1 
next  j% 
next  i% 

rmsi  *  sqr(sumi  /  cnti) 
end  sub 
•  ** 

'***  remove  the  tilt  plane  of  wavefront 
•  ** 

sub  removetilt  static 
maxp  =  -32767 
mlnp  «*  32767 
sump  =  0  :  cntp  =  0 
for  i%  =  1  to  120 
for  j%  *  1  to  120 

slope(j%,i%)  =  -32767 
next  j% 
next  1% 

for  1%  =  top%  to  bot% 
for  j%  **  xb%(l,l%)  to  xb%(2,i%) 
if  zz(j%,i%)  <>  -32767  then 

s  lope  (j%,i%)=  zz(j%,i%)  -  (coel-f  coe2  *  j  %+coe3  *  i% ) 
if  maxp  <  slope(j%,l%)  then  maxp  =  slope(j%,l%) 
if  minp  >  slope(j%,i%)  then  minp  =  slope(j%,i%) 
sump  *  sump  +  slope(  j%,  i%) ‘'2 
cntp  =  cntp  +  1 
end  if 

next  j% 
next  i% 

rmsp  sqr(sump  /  cntp) 
end  sub 
•  ** 

•***  compute  the  coefficient  of  tilt  plane 
•  ** 

sub  tiltterm  static 
n  *  0 
X  =  0 
y  -  0 
xy  »  0 
XX  «•  0 
yy  “  0 
w  *»  0 
WX  *  0 
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wy  s»  0 

for  i%  “  top%  to  bot% 
for  j%  *  xb%(l,i%)  to  xb%(2,l%) 
if  Z2(j%,i%)  <>  -32767  then 
iii  =  i% 
jjj  =  j% 

zzz  »  zz(j%,i%) 
n=n+l 
x=x+jjj 
y*y+iii 
xx=xx+j j j*2 
yy*yy+iil*2 
xy=xy+iii*j j j 
w=w+zzz 
wx=wx+zzz*j j j 
wy=wy+z z z * i i i 
end  If 
next  j% 
next  1% 

det=N*xx*yy+x*xy*y+x*y*xy-xx*y*y-N*xy*xy-x*x*yy 
coel= (w*xx*yy+wx*xy*y+x*xy*wy-xx*y*wy-w*xy*xy-wx*yy*x) /det 
coe2= (N*wx*yy+x*y*wy+w*xy*y-wx*y*y-N*wy*xy-w*x*yy ) /det 
coe3= (N*xx*wy+x*w*xy+x*y*wx-y*w*xx-N*xy*wx-x*x*wy) /det 
end  sub 


•  ** 

•***  display  the  result. 

•  ** 

sub  display  static 
call  Inltp 
call  plot(slope()) 

ys  *  2  :  xs  “  4.5  :  xsp  =  (640  -  xs  *  120) /2  :  ysp  =  50  -  top 

*  ys 

call  contour (slope 0 , lll%,rrr%,top%,bot%,xs,ys,xsp,ysp) 
key$  = 

WHILE  key$  = 

key$  *  INKEY$ 

WEND 

call  inltl 
call  plot(zz()) 

ys  “  2  :  xs  =  4.5  ;  xsp  =  (640  -  xs  *  120) /2  ;  ysp  =  50  -  top 

*  ys 

call  contour ( z z ( ) , 1 1 1 % , rrr % , top% , bot % , xs , ys , xsp , ysp ) 
key$  =  •••• 

WHILE  key$  »  "•* 
key$  -  INKEY$ 

WEND 
end  sub 


•  ** 

'***  hidden  line  remove  routine. 
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•  ** 

sub  hlddenllne(xl,x2,yl,y2,ml,n2,ininl,inm2,col%)  static 

• -  check  the  hidden  line - 

If  yl<ml  then 
If  y2<in2  then 
ininl*yl 
miii2*y2 

line  (xl,yl) -(x2,y2) , collet (col%) 
else 

in=(y2-yl)  /(X2-X1) 
b=yl-in*xl 

d=(yl-inl)  *  (x2-xl)  /  (in2-y2+yl-ml) 

x=xl+d 

y=in*x+b 

line  (xl,yl) -(x,y) ,colls%(col%) 
ininl*yl 
inin2-m2 
end  If 

elself  y2<in2  then 
in*(yl-y2)/(xl-x2) 
b*yl-m*xl 

d-(inl-yl)  *  (x2--xi)  /  (y2-m2+inl-yl) 

x*xl+d 

y«*in*x+b 

line  (x,y) -(x2,y2) ,colls%(col%) 
iiunl=inl 
inin2*y2 
else 

mmlasinl 
iiun2»in2 
end  If 
end  svYb 

•  ** 

•***  3  Dimensional  hidden  line  removed  display 
•  ** 

sub  plot(arr(2))  static 
els 

county!  ■  0 

for  1%  *  bot%  +  steprt  to  top!  step  -steprt 
prehlgh  =  0 

for  j%*lll%  to  rrr%  step  1 

If  arr(j%,l%)  <>  -32767  then 
high  =  arr ( j%, 1%) -min 
y2*hlgh*2Scale-county%*yshlft+ydc 
else 

y 2*-county %  *  y sh 1 f t+yde 
high  “  0 
end  If 

If  high  <>  0  then 

coll%  «  high  /  levellt 
coll%  «  coll%  mod  7+1 
else 
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If  prehigh  <>  0  then 
coll%  *  prehigh  /  levell% 
coll%  “  coli%  mod  7+1 
else 

coli%  =  8 
end  if 
end  if 

x2=' j  %  *xscale+xdc 
if  county%=0  then 
mp(  j0%)=*y2 

if  j%  >  111%  then  line  (xl,yl)-(x2,y2) ,colis%(coli%) 
else 

if  j%  >  111%  then 
maxl=m( j%-l) 
max2=m( j%) 

call  hlddenl ine ( xl , x2 , y 1 , y2 , maxi , max2 , maxpl , maxp2 , coli% ) 

”*P(  j^~l)*®2ixpl 

mp  ( j  % )  ^axp2 
end  if 
end  if 
xl=x2 ; yl=y2 
prehigh  =  high 
next  j% 

county%  =  county%  +  1 
next  i% 

for  j%=l  to  480 
m(j%)«mp(j%) 
next  j% 
call  drawmark 
call  boxf 
call  sec(arr()) 
end  sub 

•  ** 

•***  contour  display 

t  ** 

sub  contour  (arrd(2) ,xf%,xe%,yf%,ye%,xs,ys,xsp,ysp)  static 
els 

for  i  *  1  to  120 

x(i)  =  xsp  +  (1  -  1)  *  xs 
next  i 

for  1  »  1  to  120 

y(i)  =  ysp  +  (i  -  1)  *  ys 
next  i 

if  min  <*  0  and  max  >=  0  then 
exlevel  «  1 
else 

exlevel  **  0 
end  if 

locate  22,64 

print  levell%  chr$(248);"/  LEVEL" 
call  boxl 
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nalnp: 

for  i  =  0  to  150 

z(i)  =  min  +  i  *  levell% 
if  z(i)  >  max  then 
level  *  i 
i  =  200 
end  if 
next  i 

for  j=ye%  -  stpno%  to  yf%  step  -stpno% 

if  j<top%  or  j>=bot%  then  goto  noneinboxl 
for  i=xf%  to  xe%  -  stpno%  step  stpno% 
if  i<xb%(l,j)  or  i>xb%(2,j)  -  1  or  i  <  xb%(l, j+stpno%)  or  i  > 
xb%(2, j+stpnoi)  -  1  then  goto  noneinbox 

if  arrd(i,j)  =  -32767  or  arrd ( i+stpno% , j )  =  -32767  or 

arrd(i, j+stpno%)  =  -32767  or  arrd(i+stpno%, j+stpno%)  =  -32767  then 
goto  noneinbox 

locate  23,10 

if  (arrd(i, j) <arrd(i, j+stpno%) )  then 
dmin=arrd ( i , j ) 
else 

dmin=arrd ( i , j  +stpno% ) 
end  if 

if  arrd(i+stpno%, j)<dmin  then  dmin=arrd ( i+stpno% , j ) 
if  a r r d  (  i  +  s t p n o %  ,  j + s t p n o % )  < d m i n  then 

dmin=arrd ( i+stpno% , j  +stpno% ) 

iii  arrd(i, j) >arrd(i, j+stpno%)  then 
dmax-arrd(i, j) 
else 

dmax-arrd ( i , j  +stpno% ) 
end  if 

if  arrd ( i+stpno% , j ) >dmax  then  dmax=arrd ( i+stpno% , j ) 
if  a  r  r  d  (  i  +  s  t  p  n  o  %  ,  j -t- s  t  p  n  o  %  )  >  d  m  a  X  then 

dmax=arrd ( i+stpno% , j  +stpno% ) 

if  dmax<z(0)  or  dmin>z (level-1)  then  goto  noneinbox 
for  k=0  to  level  -  1 

if  z(k)<dmln  or  z(k)>dmax  then  goto  noneintri 
for  m=4  to  0  step  -1 
if  m>0  then 

h(m)*arrd(i+stpno%*im{m-l) , j+stpno%*jm(m-l) ) -z (k) 
xh(m)=x(i+stpno%*lm(m-l) ) 
yh(m)=y ( j+stpno%* jm(m-l) ) 
end  if 
if  m=0  then 

h(0)  =  (h(l)+h(2)+h(3)-i-h(4))/4 
xh(0)=(x(i)+x(i+stpno%) ) /2 
yh(0)=(y(j)+y(j+stpno%))/2 
end  if 

if  h(m)>0  then 
lsh(m)»2 

elseif  (h(m)<0)  then 
lsh(m)«0 
else 

ish(m)=l 
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end  if 
next  m 
for  m=l  to  4 
inl=in :  in2=0 :  in3=m+l 
if  aiS^S  then  ni3=l 

casetype-clnt(castab(lsh(nl) ,  ish(in2) ,  lsh(in3) ) ) 

If  casetype=0  then 
goto  caseO 
end  If 

on  casetype  goto 
easel , case2 , case3 , case4 , cases , case6 , case? , cases , case9 
easel: 

xl=xh(inl)  :yl=yh(inl)  :x2-xh(in2)  :y2=yh(in2) 
goto  drawlt 
case2 : 

xl-xh(in2)  :yl=yh(in2)  :x2=xh(in3)  :y2=yh(in3) 
goto  drawit 
easel : 

xl=xh(in3)  :yl=yh(in3)  :x2=xh(inl)  :y2=yh(inl) 
goto  drawlt 
case4 : 

xl=xh(inl)  :yl=yh(inl) 

x2=(h(in3)  *xh(m2)  -h(in2)  *xh(in3) )  /  (h(in3)  -h(in2) ) 
y2=(h(m3)  *yh(m2)  -h(m2)  *yh(m3) )  /  (h(m3)  -h(m2) ) 
goto  drawlt 
cases: 

xl»xh(in2)  :yl=yh(iB2) 

x2=  (h  (ml )  *xh (m3 )  -h (m3 )  *xh (ml) )  /  (h (ml)  -h (m3 ) ) 
y2=(h(ml)  *yh(m3)  -h(m3)  *yh(ml) )  /  (h(ml)  -h(m3) ) 
goto  drawlt 
cases : 

xl=xh(m3) :yl=yh(m3) 

x2=(h(m2) *xh(ml) -h(ml) *xh(m2) ) / (h(m2) -h(ml) ) 
y2=(h(m2) *yh(ml)-h(ml) *yh(m2) ) / (h(m2) -h(ml) ) 
goto  drawit 
case?: 

xl=(h(m2) *xh(ml)-h(ml) *xh(m2) ) / (h(m2) -h(ml) ) 
yl=(h(m2) *yh(ml)-h(ml) *yh(m2) ) / (h(m2) -h(ml) ) 
x2=>(h(m3)  *xh(m2)  -h(m2)  *xh(m3) )  /  (h(m3)  -h(m2) ) 
y2=(h(m3)  *yh(m2)  -h(m2)  *yh(m3) )  /  (h(m3)  -h(m2) ) 
goto  drawit 
cdS68 • 

xl=(h(m3)  *xh(ra2)-h(ro2)  *xh(m3) )  /  (h(m3) -h(m2) ) 
yl=(h(m3)  *yh(m2)-h(m2)  *yh(m3) )  /  (h(m3)  -h(m2) ) 
x2=(h(ml)  *xh(m3)  -h(m3)  *xh(rol) )  /  (h(ml)  -h(m3) ) 
y2=(h(ml)  *yh(m3)-h(m3)  *yh(ml) )  /  (h(ml)  -h(m3) ) 
goto  drawit 
case9: 

xl®(h(ml) *xh(m3)-h(m3) *xh(ml) ) / (h(ml) -h(m3) ) 
yl=(h(ml) *yh(m3) -h(m3) *yh(ml) ) / (h(ml) -h(m3) ) 
x2=(h(m2) *xh(ml)-h(ml) *xh(m2) ) / (h(m2) -h(ml) ) 
y2=(h(m2) *yh(ml) -h(ml) *yh(m2) ) / (h(ro2) -h(ml) ) 
drawit: 

xl  =  Int(xl) 
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x2  *  int(x2) 
yl  «  Int(yl) 
y2  =  int(y2) 
line  (xl,yl)-(x2,y2) 
caseO :  next  m 

nonelntri :  next  k 
nonelnbox:  next  i 
nonelnboxl:  next  j 
10202  end  sub 


•  ** 

•*** •***  draw  scale  index 
•  ** 

sub  drawmark  static 

unitx  =535  ;  unity  =45 
unitsize  =  mark!  *  zscale 

line (unitx, unity ) - (unitx, unity-unitsize) 

line (unitx-3 , unity) - (unitx+3 , unity) 

line (unitx-3 , unity-unitsize) - (unitx+3 , unity-unitsize) 

locate  4,70 

print  inark%;chr$(248)  ; 
end  sub 


•  ** 

•***  outline  box 
•  ** 

sub  boxf  static 
locate  1,8 
print  title$ 
line (0,0) -(639, 349) , ,b 
line(3, 17) -(636,346) , ,b 
llne(0,14)-(639,14) 
locate  14,50 
print  "RMS";rins  ; 
locate  14,65 
print  "P-V”;pv  ; 
end  sub 
•  ** 

•***  outline  box 
•  ** 

sub  boxl  static 
locate  1,8 
print  title$ 
line (0,0) -(639, 349) , ,b 
line (3, 17) -(636, 346) , ,b 
llne(0,14)-(639,14) 
locate  23,65 
print  "RMS";nns  ; 
locate  24,65 
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print  "P-V'';pv  ; 

end  sub 

•  ** 

•***  sectional  display  routine 
•  ** 

sub  sec(arr(2))  static 
yff  =  290 

yee  =  200 

nunff  =  0 

xff  *  111%  *  xscale  +  xdc 

xee  =»  rrr%  *  xscale  +  xdc 

mmffxc  =  (xee  -  xff)  /  (rrr%  -  111%)  *  1.190625 
degsc  =  85  /  (max  -  min  +  deginc  *  .99) 

if  top%  mod  5  then 

ydfr  =  (top%  /  5  +  1)  *  5 
else 

ydfr  =  top% 
end  if 

ydlr  =  bot%  /  5  *  5 

for  1%  =  ydfr  to  ydlr  step  20 

XX  =  (rrr%  +5)  *  xscale  +  xdc 
yy  =>  (i%  -  bot%)  /  stepr%  *  yshift  +  ydc 
line(xx,yy) -(xx+15,yy) 
row%  «  yy  /14  +  1 
col%  =  XX  /  8  +  3 
locate  row%,col% 
print  i%  ; 
next  1% 
secno%  =  0 

do  while  secno%  <>  999 
secno%  =  0 

do  until  secno%  =  999  or  secno%  >=  top%  and  secno%  <=  bot% 
locate  23,5 

print  "Key  in  trace  no."?top%;"-";bot%;  "  999  to  quit  " 

•I  . 

I 

locate  23,50 
input  secno% 

if  secno%  <>  999  and  (secno%  <  top%  or  secno%  >  bot%)  then 
print  chr$(7)  ; 
loop 

if  secno%  <>  999  then 
call  clearsec 
locate  23,35 
print  "TRACE  #";secno% 
locate  22,6 
print  "(mm)"  ; 

for  degff=0  to  (max-mln)  +  deginc  *.99  step  deginc 
yy  =  -degff  *  degsc  +  yff 
line(xff ,yy)-(xee,yy) , , ,&hcccc 
locate  yy  /  14  +  1,68 
print  degff ;chr$ (248)  ; 
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next  degff 

for  mmff  =  0  to  99  step  10 
XX  =  mmff  *  mmffxc  +  xff 
line(xx,yff ) -(xx,yy) , , ,&hcccc 
nunia%  =  minff 

if  minint  mod  20  =  0  then 
locate  yff/14+l,xx/8 
print  mmm%  ; 
end  if 
next  mmff 

line(xee,yff)-(xee,yy)  ,,,&hcccc 

first  =  1 
sum  =  cnt  =  0 

peak  =  -32767  :  valley  =  32767 

for  j%  *  111%  to  rrr% 

if  arr( j%,secno%)  <>  -32767  then 
XX  =  j%  *  xscale  +  xdc 
yy  =  (min-arr ( j%,secno%) )  *  degsc  +  yff 
if  peak  <  arr ( j%,secno%)  then  peak  =  arr ( j%,secno%) 
if  valley  >  arr ( j%,secno%)  then  valley  =  arr ( j%,secno%) 
sum  =  sum  +  arr  ( j%,secno%)  “'2 
cnt  =  cnt  +  1 

if  first  <>  1  then  line(xx,yy) -(xxo,yyo) , 14 
first  =  0 

XXO  =  XX 

yyo  =  yy 
end  if 
next  j% 
locate  23,50 

print  ''RMS'';sqr  (sum/cnt)  ; 
locate  23,65 

print  "P-V”;peak-valley  ; 
key$  =  '•« 

WHILE  key$  = 

key$  =  INKEY$ 

WEND 

end  if 
loop 
end  sub 


I  ** 

•***  clean  clipregin 
<  ** 

sub  clearsec  static 
locate  15,2 
print  ” 

ti  . 

$ 

locate  16,2 
print  " 

It  . 

0 


locate  17,2 
print  " 


ft 


9 
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locate  18,2 
print  " 

H 

locate  19,2 
print  " 

II 

locate  20,2 
print  " 

H 

locate  21,2 
print  '• 

II 

locate  22,2 
print  " 

II 

locate  23,2 
print  ” 

II 

locate  24,2 
print  " 

II 

end  sub 


I*  ************************* 

*  ACQ.C  * 

************************* 

*/ 

/*  Get  image  by  using  EPIX  and  MDC3-BMCC  */ 

#include  <process.h> 

#include  <string.h> 

#include  <stdio.h> 

#include  <stdlib.h> 
jfinclude  <signal.h> 

#include  <dos.h> 

#include  "svobj.h”  /*  function  prototypes  */ 

#include  "pxmp.h” 

int  xdiin,ydiin,  images  ; 

long  int  buf, speed  ; 

char  one[80] ,tot[80] , tmp[80] ,spd[80]  ; 

ma in ( msgn , msgc ) 
int  msgn  ; 
char  *rosgc [ ]  ; 

{  long  int  buf s, buf e  ; 

int  i, ii,xwid,yintl,j ,k, step, steps  ; 


if  (msgn  <  2) 

{  printf ("\n  'Getiraage  XDIM  YDIM  BUFNO  IMAGES  STEPS  SPEED'  to 
run  " )  ; 

printf  ("\nFunct ion  :  Get  image  by  using  EPIX  and  MDC3-BMCC") 

9 
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printf  (''\nParaineter!  ")  ; 

printf(”\n  Bufno  :  Starting  buffer  no  of  images  in  SV 
board”)  ; 

printf (”\n  FILENM  :  IMAGE  NAME  to  be  save”)  ; 
exit (99)  ; 

} 

xdim  =  ydim  =  240  ; 

buf  =»  atol  (msgc[l] )  ; 

step  =  63  ; 

steps  =  step  *  5  ; 

images  =  5  ; 

ii  =s  images  -  1  ; 

sprintf (one, ”+T%d,S”, step)  ; 

sprintf (tot, ”-T%d,S”, steps)  ; 


opensi(240,3,0,240,-l,0)  ; 
openoriel()  ; 
bufs  =  buf  ; 
bufe  =  bufs  +  4  ; 

for  (i  =  0  ;  i  <  5  ;  i++,buf++) 

{ 

tocoml(one)  ; 

for  (j  =  -32767  ;  j  <  32767  ;  j++) 
for  (k  =  0  ;  k  <  30  ;  k++)  ; 
pxd_video ( ' z ' , buf )  ; 
pxd_video ( ' p ' , buf )  ; 

} 

tocoml(tot)  ; 

/*  pxd  video( 's' ,0L)  ; 

i  f 

( ( er rno=pxio8_f mbwr i te ( pxd_def image (OL, 0,0, -1,-1) , bufs, bufe, msgc [ 
2],0L,1,1))  !=  0) 

printf  (''\nError;  %s'',pxerrnomesg(errno) )  ;  */ 

pxd_close ( )  ; 

} 


opensi (xdim, xwid,xoff, ydim, yoff, intla) 
int  xdim,xwid,xoff , ydim, yoff , intla  ; 

{  do_signal()  ; 

if  (pxd  svopen(0xD,  0x300,  5,  0xD)<0) 

{  prTntfC'Open  Silicon  Video  Error  \n”)  ; 
exit(99)  ; 

> 

pxd_vidgenlock(l, 1)  ; 

pxd_vidparm(xdiro,xwid,xoff , ydim, 1, yoff , intla)  ; 

} 

sigintO 

{ 

int  abort; 

printf  (''\nlnterrupt  -  closing  SILICON  VIDEO\n”) ; 
pxd_close() ; 
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exit(l) ; 

} 

sigfpeO 

{ 

int  abort; 

printf  ("\nFloating  point  exception  -  closing  SILICON  VIDEO\n*') ; 
pxd_close ( ) ; 

/*  setvideoinode(_TEXTC80)  ;  */ 

exTt(l) ; 


do_signal() 

{ 

signal (SIGINT,  sigint) ; 
signal (SIGFPE,  sigfpe) ; 


openorielO 
{  int  i ; 

union  REGS  in, out  ; 

char  lnit[20]  =  "\x003\x00aA\O”  ; 

(in) .x.dx  =  0  ; 

(in).h.al  =  0xe2  ; 

(in) .h.ah  =  0  ; 

Int86(0xl4,&in,&out)  ; 
tocomlC'AEH”)  ; 

} 

tocoml (string) 
char  *string; 

{ 

int  j , ch ; 

for  (  ;  *string  1=  '\0'  ;  string++) 

{  for  (j=-20000;  j  <  2000;  j++) ;  /*  the  maxliaum  speed  of 

1200  */ 

ch  =  *string; 
outp  (0x3f8,ch); 

} 

> 
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I  It  1l  11  h  h  h  h  h  It  h  It  h  IHt  It  h  it  It  h  h  h  h  Ic  h  h 

•  *  DEMO. BAS  * 

I  AA  AAA  AA  A* 


'This  is  the  complete  demo  package  for  the  multitrace  system 


COMMON  SHARED  selkey  'declare  command  variable 
const 
thold%=2 00 ,  cutr %=10 ,  cutc%=2 0 ,  ton%®l ,  tof  f  %=0 ,  center%=12  0 ,  size%=:2 4  0 
, initord%=0 

const  scalepi=. 5729583, smth%=l 
const  shift%  =  3 

const  stepr%«*3,  xscale=4.2,yscale*2.5,yshift=3.8,stpno%  =  3 

const  xdc  =  20,ydc  =  170,111%  =  l,rrr%  =  120 

common  shared  top%,bot%,ib%,rb% 

common  shared  max,min,bu£% 

common  shared  row%,set% 

common  shared  coel,coe2,coe3 

common  shared  maxp,minp, maxi, mini, rmsp,rmsi,pvp,pvi,rms,pv 
common  shared  levellt, zscale,deginc,mark%, title$ 
rem  $dynamic 

dim  shared 

Obj (240) ,C%(240) ,r%(240) ,ZZZ(120,5) ,tz(360) ,ttz(360) , shif tpts% (5) 

dim  shared 

order%(240) ,slope(120, 120) , zz (120, 120) ,xb% (2, 120) ,bb%(2,5) , jff%(5 

) , jee%(5) ,filenm$(5) 

dim  shared  colis%(8) 

dim  shared  m(480) ,mp(480) 

dim  shared 

z(150) ,h(9) ,ish(9) ,xh(5) ,yh(5) ,im(9) ,x(120) ,y(120) , jm(5) ,castab(3 
,3,3) 

rem  $statlc 

data 


0,0, 8, 0,0, 2, 5, 1,7, 6, 9, 2, 3, 4, 5, 6, 0,3, 4, 7, 1,0, 1,8, 4, 3, 0,9, 0,1, 2, 3, 9 


screen  9 


call  printmenu 
CALL  initial 

DO  UNTIL  selkey  =  27  'ESC  key  to  exit 

program 

CALL  kbhit  'wait  for  key  hit  from 

operator 

SELECT  CASE  selkey 

CASE  59  'have  the  stage  moved 

to  each 

CALL  send(''MR+l8750'')  'aperture  for  testing 

scanning 

call  stopped 

CALL  send(''MR+18750'') 

call  stopped 

CALL  send(''MR+18750") 
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•  ************************* 
•  *  DEMO. BAS  * 

I  ************************* 


'This  Is  the  complete  demo  package  for  the  multitrace  system 


COMMON  SHARED  selkey  'declare  command  variable 
const 
thold%=2  00 , cutr%=10 , cutc%=20 , ton%=l , tof f %=0 , center%*120 , slze%=24  0 
, lnltord%=0 

const  scalepl=. 5729583, smth%=l 
const  shlft%  -  3 

const  stepr%=3,  xscale=4.2,yscale=2.5,yshlft=3.8,stpno%  -  3 

const  xdc  »  20,ydc  *  170,111%  =  l,rrr%  =  120 

common  shared  top% , bot% , lb% , rb% 

common  shared  max,mln,buf% 

common  shared  row%,eet% 

common  shared  coel,coe2,coe3 

common  shared  maxp,mlnp, maxi, mini, rmsp,rmsl,pvp,pvl,rms,pv 
common  shared  levellt, zscale,deglnc,mark%, tltle$ 
rem  $dynamlc 

dim  shared 

obj (240) ,c%(240) ,r%(240) ,zzz(120,5) ,tz(360) ,ttz(360) , shlftptst (5) 

dim  shared 

order%(240) ,slope(120, 120) ,zz(120,120) ,xb%(2,120) ,bb%(2,5) , jff%(5 

) , jee%(5) ,fllenm$(5) 

dim  shared  colls%(8) 

dim  shared  m(480) ,mp(480) 

dim  shared 

z (150) ,h(9) , lsh(9) ,xh(5) ,yh(5) , lm(9) ,x(120) ,y (120) , jm(5) ,castab(3 
,3,3) 

rem  $statlc 

data 


,6,7,- 


screen  9 


call  prlntmenu 
CALL  Initial 
DO  UNTIL  selkey  *  27 
program 

CALL  kbhlt 
operator 

SELECT  CASE  selkey 
CASE  59 
to  each 

CALL  send("MR+18750") 
scanning 

call  stopped 

CALL  send("MR+18750'') 

call  stopped 

CALL  send("MR+18750'') 


•ESC  key  to  exit 
’wait  for  key  hit  from 


•have  the  stage  moved 
'aperture  for  testing 
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call  stopped 
CALL  send('’MR+18750") 
call  stopped 
CALL  send  (••MR-75000") 
call  stopped 
call  prlntmenu 
CASE  60 
CLS 

LOCATE  12,10 

print  "  ACQUISITION  SUBAPERTURE  #1 

n 

shell ("acq  5")  •call  acg.exe  to  get  the 

image 

call  waiting 
CALL  send("MR+18750") 
call  stopped 
LOCATE  12,10 

print  "  ACQUISITION  SUBAPERTURE  §2 


shell ("acq  10") 
call  waiting 
CALL  send("MR+18750") 
call  stopped 
LOCATE  12,10 

print  "  ACQUISITION  SUBAPERTURE  #3 

n 

shell ("acq  15") 
call  waiting 
CALL  send("MR+18750") 
call  stopped 
LOCATE  12,10 

print  "  ACQUISITION  SUBAPERTURE  #4 

n 

Shell ("acq  20") 
call  waiting 
CALL  send('^MR+18750") 
call  stopped 
LOCATE  12,10 

print  "  ACQUISITION  SUBAPERTURE  #5 

n 

shell ("acq  25") 
call  waiting 
CALL  send ("MR-75000") 
call  stopped 
call  main 
call  prlntmenu 
CASE  ELSE 
END  SELECT 
LOOP 

els 

END 


»***  Dummy  waiting  subroutine. 
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•  ** 

sub  waiting  static 

for  ii  =  -32767  to  32767 
for  jj  “  0  to  1 
next  jj 
next  ii 
end  sub 


t  ** 

•***  main  process,  contains  computation  ,  and  display. 

•  ** 

sub  main  static 
call  init 
els 

locate  12,1 

print  ••  COMPUTING  ... 

H 

call  waiting 
call  display 
end  sub 


t  ** 

•  ***  wait  for  asterisk  to  make  sure  command  is  accept  by  controller 

•  ** 

SUB  ASTERISK  STATIC 

DO  UNTIL  INPUT$(1,  #6)  = 

LOOP 

END  SUB 

•  ** 

'***  initial  process 

•  ** 

SUB  initial  STATIC 

OPEN  '*COM2:9600,N,8,1,RS,CS,DS,CD''  FOR  RANDOM  AS  #6  'OPEN  SERIAL 
PORT 

CALL  send ("MED”)  'disable  mnemonics 

CALL  send("MFE")  'enable  move  finish 

CALL  send("VI600")  'set  initial  velocity 

CALL  send("VF6000")  'set  final  velocity 

CALL  8end("AC1001")  'set  acceleration  and 

deceleration 

END  SUB 

•  ** 

•***  accept  key  hit  from  keyboard 
•  ** 

SUB  kbhit  STATIC 
key$  - 
DO 

WHILE  key$  =  "" 
key$  »  INKEY$ 

WEND 
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IP  LEN(key$)  >  1  THEN 

skey$  =  MID$(key$,  2,  1) 

ELSE 

skey$  *  MID$(key$,  1,  1) 

END  IF 

selkey  »  ASC(skey$) 

LOOP  UNTIL  key$  <>  ”” 

END  SUB 

•  ** 

•***  subroutine  print  the  operation  menu  onto  screen 

•  ** 

SUB  printmenu  STATIC 
CLS 

LOCATE  3,1 

print  •*  CYLINDRICAL  MIRROR  MULTITRACE  MEASUREMENT 

H 

print 
print 
print  "" 
print 

print  "  FI  :  TEST  MIRROR  SCANNING” 

print 

print  •*  F2  :  ACQUISITION  AND  MEASUREMENT” 

print  ”” 
print  ”” 

print  ”  ESC:  EXIT” 

END  SUB 

•  ** 

»***  subroutine  for  sending  command  to  controller 

I  ** 

SUB  send  (a$)  STATIC 

PRINT  #6,  a$  'send  command  to 

controller 

CALL  ASTERISK  'wait  until  asterisk 

is  return 
END  SUB 

•  ** 

'***  wait  until  process  finish 

I  ** 

SUB  STOPPED  STATIC 

DO  UNTIL  INPUT$(1,  #6)  -  ”F" 

LOOP 

for  ii  *  -32767  to  32767 
for  jj  *  0  to  1 
next  jj 
next  ii 

END  SUB 


•  ** 

'***  initialize  parameter. 
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•  ** 

sub  init  static 
colls%(l)  *  10 
colls%(2)  =  3 

colls%(3)  a=  9 
colls%(4)  **  1 

colls%(5)  =  4 

colis%(6)  s  12 
colis%(7)  *  14 
colls%(8)  =  15 

iin(0)*'0  ;  iin(l)=l  :  itn(2)=l  :  im(3)“0 
jm(0)*0  :  jm(l)»=0  :  jin(2)=l  :  jin(3)=l 
for  k“0  to  3 
for  j*0  to  3 

for  1*0  to  3 
read  castab(k, j , 1) 
next  1 
next  j 
next  k 

open  "deinor.dat"  for  Input  as  #1 
open  "demol.dat"  for  Input  as  #2 
Input  #1, top%,bot%,inaxp,inlnp 
Input  #2,topt,bot%,inaxl,inlnl 
for  1%  *  1  to  120 

for  j%  =  1  to  120 
zz(j%,l%)  =■  -32767 
slope(j%,l%)  =  -32767 
next  j% 
next  It 

sump  =  0  :  suml  =  0  :  cntp  =  0  :  cntl  =  0 
for  1%  =  topt  to  bott 

Input  #l,dmy%,xb%(l,l%) ,xb%(2,l%} 
Input  #2,dmy%,xb%(l,l%) ,xb%(2,l%) 
for  j%  *  xb%(l,l%)  to  xb%(2,l%) 

Input  #1, slope ( j%, 1%) 

Input  #2,zz(j%,l%) 

If  slope(j%,l%)  <>  -32767  then 
sump  =  sump  +  slope(  j%,  1%)  "^2 
cntp  =  cntp  +  1 
end  If 

If  zz(j%,l%)  <>  -32767  then 
sum!  =  suml  +  zz(j%,l%)'‘2 
cntl  -  cntl  +  1 
end  If 
next  j% 
next  1% 

rmsp  *  sqr(sump  /  cntp) 
rmsl  *  sqr(suml  /  cntl) 
close 
end  sub 


•***  Initialize  parameter  of  slope  display. 
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sub  Inltp  static 

title$="CYLINDRICAL  MIRROR  :  MULTITRACE  MEASUREMENT  /  SLOPE 
FUNCTION" 

levell%  =  12 
zscale  =  -.3 
deginc  *  36 
mark!  -  60 
for  i%  =  1  to  480 
m(i%)  =  32767 
inp(i%)  =  32767 
next  1% 
max  s  maxp 
min  =  minp 
rms  s  rmsp 
pv  =  maxp  -  mlnp 
end  sub 

•  ** 

»***  initialize  parameter  of  wavefront  display. 

•  ** 

sub  Inltl  static 

title$="CYLINDRICAL  MIRROR  :  MULTITRACE  MEASUREMENT  /  WAVEFRONT 
FUNCTION" 

levell%  =  180 
zscale  =»  -.005 
deginc  -720 
mark%  =  1800 
for  i%  =  1  to  480 
m(l%)  =  32767 
mp(i%)  =  32767 
next  i% 
max  a  maxi 
min  a  mini 
rms  a  rmsi 
pv  a  maxi  -  mini 
end  sub 


•  ** 

'***  display  the  result. 

» ** 

sub  display  static 
call  initp 
call  plot(slope() ) 

ys  a  2  :  xs  a  4,5  :  xsp  =  (640  -  xs  *  120) /2  :  ysp  a  50  -  top 

*  ys 

call  contour (slope 0 , 111%, rrr%, topi, bot%,xs,ys, xsp, ysp) 
key$  s  nn 

WHILE  key$  * 

key$  a  INKEY$ 

WEND 

call  inltl 
call  plot (zz  0 ) 
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ys  *  2  :  xs  »  4.5  :  xsp  *  (640  -  xs  *  120) /2  :  ysp  =  50  -  top 

*  ys 

call  contour (zz() , lll%,rrr%,top%,bot%,xs,ys,xsp,ysp) 
key$  * 

WHILE  key$  = 

key$  =  INKEY$ 


WEND 
end  sub 


•  ** 

’***  hidden  line  remove  routine. 

•  ** 

sub  hlddenllne(xl,x2,yl,y2,ml,m2,mml,inm2,col%)  static 

* -  check  the  hidden  line - 

If  yl<ml  then 
if  y2<m2  then 
mml=yl 
mm2=y2 

line  (xl,yl) -(x2,y2) ,colis% (col%) 
else 

m=(y2-yl)/(x2-xl) 

b=yl-m*xl 

d=(yi-roi) * (x2-xl) / (m2-y2+yl-ml) 

x=xl+d 

y=m*x+b 

line  (xl,yl)-(x,y) ,colis%(col%) 
ramls=yl 
mm2=m2 
end  If 

elself  y2<m2  then 
m=(yl-y2)/(xl-x2) 
b=yl-m*xl 

d='(ral-yl)  *  (x2-xl)  /  (y2-m2+ml-yl) 

x=xl+d 

y=m*x+b 

line  (x,y) -(x2,y2) ,colis%(col%) 
mml=ml 
mm2=y2 
else 

mml=ml 
inm2=m2 
end  if 

end  sub 

•  ** 

I***  3  Dimensional  hidden  line  removed  display 

•  ** 

sub  plot(arr(2))  static 
els 

county!  *=  0 

for  i%  =  bot%  +  stepr%  to  top!  step  -stepr! 
prehigh  =  0 
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for  j%=lll%  to  rrr%  step  1 

if  arr(j%,l%)  <>  -32767  then 
high  *  arr ( j%, i%) -min 
y2=high*zscale-county%*yshift+ydc 
else 

y2*-county%*yshift+ydc 
high  *  0 
end  if 

if  high  <>  0  then 

coli%  =  high  /  levell% 
coli%  *  coli%  mod  7+1 
else 

if  prehigh  <>  0  then 
coli%  =  prehigh  /  levellt 
coli%  =  coli%  mod  7+1 
else 

coll%  =  8 
end  if 
end  if 

x2=j%*xscale+xdc 
if  county %=0  then 
mp( j0%)=y2 

if  j%  >  111%  then  line  (xl,yl)-(x2,y2) ,colis%(coli%) 
else 

if  j%  >  111%  then 
maxl*m( j%-l) 
max2=m(j%) 

call  hiddenline ( xl , x2 , yl , y2 , maxi , max2 , maxpl , maxp2 , col i% ) 
mp( j%-l)=maxpl 
mp( j%)*maxp2 
end  if 
end  if 
xl=x2 :yl=y2 
prehigh  =  high 
next  j% 

county%  =  county%  +  1 
next  i% 

for  j%=l  to  480 
m(j%)=mp( j%) 
next  j% 
call  drawmark 
call  boxf 
call  sec(arr()) 
end  sub 

•  ** 

•***  contour  display 
•  ** 

sub  contour  (arrd(2) ,xf%,xe%,yf%,ye%,xs,ys,xsp,ysp)  static 
els 

for  i  »  1  to  120 

x(l)  *  xsp  +  (1-1)  *  xs 
next  i 

for  i  =•  1  to  120 
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y(i)  =  ysp  +  (i  -  1)  *  ys 
next  1 

if  min  <=  0  and  max  >=  0  then 
exlevel  *  1 
else 

exlevel  »  0 
end  if 

locate  22,64 

print  levell%  chr$(248);”/  LEVEL” 
call  boxl 

mainp: 

for  i  =  0  to  150 

z(i)  ®  min  +  i  *  levell% 
if  z(i)  >  max  then 
level  =  i 
i  =  200 
end  if 
next  i 

for  j=ye%  -  stpno%  to  yf%  step  -stpno% 

if  j<top%  or  j>*=bot%  then  goto  noneinboxl 
for  i=*xf%  to  xe%  -  stpno%  step  stpno% 
if  i<xb%(l,j)  or  i>xb%(2,j)  -  1  or  i  <  xb% (1, j+stpno%)  or  i  > 
xb%(2, j+stpno%)  -  1  then  goto  noneinbox 

if  arrd(i,j)  =  -32767  or  arrd(i+stpno%, j)  =  -32767  or 

arrd(i, j+stpno%)  =*  -32767  or  arrd(i+stpno%, j+stpno%)  =  -32767  then 
goto  noneinbox 

locate  23,10 

if  (arrd(i, j)<arrd(i, j+stpno%) )  then 
dmin=arrd(i, j) 
else 

dmin=arrd ( i , j  +stpno% ) 
end  if 

if  arrd(i+stpno%, j) <dmin  then  dmin=arrd ( i+stpno% , j ) 
if  a r r d  (  i  +  s t p n o %  ,  j + s t p n o % )  < d m i n  then 

droin=arrd ( i+stpno% , j  +stpno% ) 

i f  ar rd ( i , j ) >arr d ( i , j  +stpno% )  then 
dmax=arrd ( i , j ) 
else 

dmax=ar rd ( i , j  +s tpno% ) 
end  if 

if  arrd(i+stpno%, j) >dmax  then  dmax=arrd ( i+stpno% , j ) 
if  arrd(i+stpno%,j+stpno%)>dmax  then 

dmax=arrd ( i+stpno% , j+stpno% ) 

if  dmax<z(0)  or  droin>z (level-1)  then  goto  noneinbox 
for  k=0  to  level  -  1 

if  z(k)<dmln  or  z(k)>dmax  then  goto  noneintri 
for  m=4  to  0  step  -1 
if  m>0  then 

h(m)=arrd(l+stpno%*lm(m-l) , j+stpno%*jm(m-l) ) -z (k) 
xh(m)=x(i+stpno%*lm(m-l) ) 
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yh(in)=y  ( j+stpno%*jin(ro-l) ) 
end  if 
if  m=0  then 

h(0)=(h(l)+h(2)+h(3)+h(4))/4 
xh(0)=(x(i)+x(i+stpno%) ) /2 
yh(0)*(y(j)+y(j+stpno%) )/2 
end  if 

if  h(in)>''  then 
ish(in)=2 

elseif  (h(]n)<0)  then 
ish(in)=0 
else 

ish(in)=l 
end  if 
next  m 
for  in=l  to  4 
inl=in :  in2=0 :  in3=m+l 
if  in3=5  then  in3=l 

casetype=cint (castab(ish(inl) ,  ish(in2) ,  ish(in3) ) ) 
if  casetype=0  then 
goto  caseO 
end  if 

on  casetype  goto 
easel , case2 , case3 , case4 , cases , case6 , case? , cases , case9 
easel: 

xl=xh(inl)  :yl*yh(inl)  :x2=xh(in2)  :y2=yh(in2) 
goto  drawit 
case2 : 

xl=xh(in2)  :yl*=yh(in2)  :x2=xh(in3)  :y2=yh(in3) 
goto  drawit 
case3 : 

xl=xh(in3)  :yl=yh(in3)  ;x2=xh(inl)  :y2=yh(ml) 
goto  drawit 
case4: 

xl=xh(inl)  :yl=yh(inl) 

x2=(h(in3)*xh(in2)-h(in2)*xh(in3) )  /  (h(in3) -h(in2) ) 
y2=(h(Tn3)  *yh(in2)-h(in2)  *yh(in3) )  /  (h(in3)  -h(in2) ) 
goto  drawit 
cases: 

xl=xh(in2)  :yl=yh(in2) 

x2=(h(Tnl)  *xh(in3)-h(in3)  *xh(inl) )  /  (h(inl) -h(Tn3) ) 
y2=(h(inl)  *yh(in3)-h(in3)  *yh(Tnl) )  /  (h(inl)  -h(in3) ) 
goto  drawit 
cases: 

xl=xh(in3)  :yl=yh(ro3) 

x2=(h(in2)  *xh(inl)  -h(inl)  *xh(in2) )  /  (h(in2)  -h(Tnl) ) 
y2=(h(Tn2)  *yh(ml)~h(inl)  *yh(in2) )  /  {h(in2)  -h(rol) ) 
goto  drawit 
case? : 

xl=(h(in2)  *xh(Tnl)-h(inl)  *xh(in2) )  /  (h(m2)  -h(ml) ) 
yl=(h(in2)  *yh(ml)  -h(inl)  *yh(in2) )  /  (h(in2)  -h(inl) ) 
x2=(h(in3)  *xh(in2)-h(in2)  *xh(in3) )  /  (h(in3) -•h(in2) ) 
y2=  (h  (in3 )  *yh  (m2 )  -h  (ro2 )  *yh  (ro3 ) )  /  (h  (m3 )  -h  (m2 ) ) 
goto  drawit 
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cases : 

xl=(h(in3)  *xh(in2)  -h(in2)  *xh(in3) )  /  (h(in3)  -h(m2) ) 
yl=(h(in3)  *yh(m2)  -h(in2)  *yh(m3) )  /  (h(in3)  -h(m2) ) 
x2=(h(inl)  *xh(in3)  “h(in3)  *xh(Tnl) )  /  (h(inl)  -h(in3) ) 
y2=(h(ml)  *yh(in3)  -h(in3)  *yh(iol) )  /  (h(ml)  -h(in3) ) 
goto  drawit 
case9 : 

xl={h(inl)  *xh(in3)  -h(in3)  *xh(inl) )  /  (h(inl)  -h(in3) ) 
yl=(h(Tnl)  *yh(Tn3)  -h(in3)  *yh(inl) )  /  (h(Tnl)  -h(m3) ) 
x2=(h(m2)  *xh(inl)  -h(inl)  *xh(in2) )  /  (h(m2)  -h(inl) ) 
y2=(h(m2)  *yh(inl)  -h(inl)  *yh(in2) )  /  (h(in2)  -h(inl) ) 
drawit: 

xl  =  int(xl) 
x2  =  int(x2) 
yl  =  int(yl) 
y2  =  lnt(y2) 
line  (xl,yl)-(x2,y2) 
caseO:  next  m 


noneintri:  next  k 
noneinbox:  next  i 
noneinboxl;  next  j 
10202  end  sub 


•  ** 

'***  draw  scale  index 
•  ** 

sub  drawmark  static 

unitx  =  535  :  unity  =45 
unitsize  =  inark%  *  zscale 

line (unitx, unity) - (unitx, unity-unitsize) 
line (unitx-3 , unity) - (unitx+3 ,unity) 

line (unitx-3 , unity-unitsize) - (unitx+3 , unity-unitsize) 
locate  4,70 

print  inark%;chr$  (248)  ; 
end  sub 


•  ** 

•***  outline  box 
•  ** 

sub  boxf  static 
locate  1,8 
print  title$ 
line (0,0) -(639, 349)  ,  ,b 
line(3,17)-(636,346) ,  ,b 
line(0,14)-(639,14) 
locate  14,50 
print  '*RMS”;rins  ; 
locate  14,65 
print  ''P-V'’;pv  ; 
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end  sub 
•  ** 

'***  outline  box 
•  ** 

sub  boxl  static 
locate  1,8 
print  title$ 
line (0,0) -(639, 349)  ,  ,b 
line (3, 17) -(636, 346) ,  ,b 
line(0,14)-(639,14) 
locate  23,65 
print  ''RMS";rias  ; 
locate  24,65 
print  ”P-V”;pv  ; 
end  sub 

•  ** 

•***  sectional  display  routine 
•  ** 

sub  sec(arr(2))  static 
yff  =  290 

yee  =  200 

mmff  =  0 

xff  =  111%  *  xscale  +  xdc 

xee  =  rrr%  *  xscale  +  xdc 

mmffxc  =  (xee  -  xff)  /  (rrr%  -  111%)  *  1.190625 
degsc  -  85  /  (max  -  min  +  deginc  *  .99) 

if  top%  mod  5  then 

ydfr  =  (top%  /  5  +  1)  *  5 
else 

ydfr  =  top% 
end  if 

ydlr  =  bot%  /  5  *  5 

for  i%  =  ydfr  to  ydlr  step  20 

XX  =  (rrr%  +  5)  *  xscale  +  xdc 
yy  =  (i%  -  bot%)  /  stepr%  *  yshift  +  ydc 
line(xx,yy) -(xx+15,yy) 
row%  =  yy  /14  +  1 
col%  =  XX  /  8  +  3 
locate  row%,col% 
print  i%  ; 
next  i% 
secno%  =  0 

do  while  secno%  <>  999 
secno%  =  0 

do  until  secno%  =  999  or  secno%  >=  top%  and  secno%  <=  bot% 
locate  23,5 

print  "Key  in  trace  no.";top%;"-";bot%;  "  999  to  quit  " 

n  . 

I 

locate  23,50 
input  secno% 

if  secno%  <>  999  and  (secno%  <  top%  or  secno%  >  bot%)  then 
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print  chr$(7)  ; 
loop 

If  secno%  <>  999  then 
call  clearsec 
locate  23,35 
print  "TRACE  #";secno% 
locate  22,6 
print  " (mm) "  ; 

for  degff=0  to  (max-min)  +  deginc  *.99  step  deginc 
yy  =  -degff  *  degsc  +  yff 
line(xff ,yy) -(xee,yy) , , ,6hcccc 
locate  yy  /  14  +  1,68 
print  degff ;chr$(248)  ; 
next  degff 

for  mmff  =  0  to  99  step  10 
XX  =  mmff  *  mmffxc  +  xff 
line(xx,yff ) -(xx,yy) , , ,&hcccc 
mmm%  =  mmff 

if  mmm%  mod  20  =  0  then 
loc'+'e  yff /14+l,xx/8 
print  mmm%  ; 
end  if 
next  mmff 

line(xee,yff) -(xee,yy)  ,,,!thcccc 

first  =  1 
sum  =  cnt  =  0 

peak  *  -32767  ;  valley  =  32767 

for  j%  =  111%  to  rrr% 

if  arr( j%,secno%)  <>  -32767  then 
XX  =»  j%  *  xscale  +  xdc 
yy  =  (min-arr ( j%,secno%) )  *  degsc  +  yff 
if  peak  <  arr ( j%,secno%)  then  peak  =  arr ( j%,secno%) 
if  valley  >  arr ( j%,secno%)  then  valley  =  arr( j%,secno%) 
sum  =  sum  +  arr  ( j%,secno%)  ^^2 
cnt  =  cnt  +  1 

if  first  <>  1  then  line(xx,yy)-(xxo,yyo) ,14 
first  =  0 
XXO  =  XX 
yyo  =  yy 
end  if 
next  j% 
locate  23,50 

print  "RMS";sgr(sum/cnt)  ; 
locate  23,65 

print  "P-V";peak-valley  ; 
key$  =  "" 

WHILE  key$  =  "" 

key$  =  INKEY$ 

WEND 

end  if 
loop 
end  sub 
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I  ** 

•***  clean  clipregin 

•  ** 

sub  clearsec  static 
locate  15,2 
print  " 

It  • 

9 

locate  16,2 
print  ” 

It  • 

9 

locate  17,2 
print  « 

II  • 

9 

locate  18,2 
print  ” 

It  • 

9 

locate  19,2 
print  ” 

II  • 

9 

locate  20,2 
print  ” 

II  « 

9 

locate  21,2 
print  •• 

II  • 
t 

locate  22,2 
print  ” 

II  • 

9 

locate  23,2 
print  ” 

II  • 

9 

locate  24,2 
print  ” 

II  • 
t 

end  sub 
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6.0  TESTING  AND  EXPERIMENTAL  RESDLTfl 

€.1  The  Teat  Mirror 

A  cylindrical  convex  surface  was  custom  fabricated  for  testing 
purposes.  The  surface  is  characterized  as  follows: 

Dimension:  3”  wide  x  9"  long 

Radius  of  curvature:  300mm 

Surface  errors:  figure  and  macroroughness. 

The  mirror  is  mounted  on  a  tilt/rotation  manual  stage  for  alignment 
purposes.  This  in  turn  is  mounted  on  the  NEAT  scanning  stage  for 
subaperture  testing. 


€.2  subaperture  Interferoarama 

In  the  following,  we  present  a  number  of  subaperture  inter ferograms 
for  various  parts  of  the  mirror  and  various  interferometer 
alignment  conditions. 

Three  subapertures  were  chosen  along  the  mirror  surface: 

-  Subaperture  1:  Medium  surface  quality 

“  Subaperture  2:  good  surface  quality 

Subaperture  3:  Poor  surface  quality. 

For  the  subapertures,  interferograms  are  captured  with  various 
interferometer  alignment  conditions  as  follows: 

Grating  in  focus  and  mirror  aligned 

Grating  in  focus  and  mirror  angularly  misaligned  by  various 
amounts 

Grating  out-of  focus  and  mirror  aligned 

Grating  out-of  focus  and  mirror  angularly  misaligned  by 
various  amounts. 

Photographs  of  the  various  interferograms  are  shown  in  Figure  8 
with  self-explanatory  captions. 
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6.3  Subaperture  Testing:  detailed  Operation 


In  the  following  we  present  the  results  of  the  testing  program 
showing  the  details  of  the  subaperture  testing  approach. 

The  test  makes  use  of  5  subapertures.  Each  subaprture  is  2”  long 
and  the  overlap  area  between  subapertures  is  1”,  so  that  the  total 
scan  length  is  about  6”  long.  For  each  subaperture  the  phase 
function,  i.e.  the  wavefront  slope  function,  is  first  independently 
computed.  Then,  simple  first  order  least  square  fitting  is  carried 
out  in  the  overlap  area  to  create  a  continuous  function  along  the 
mirror. 

The  results  are  given  in  Figures  9  to  12: 

Figure  9  shows  photographs  of  the  5  subaperture 
inter ferograms,  separate  and  side  by  side 

Figure  10  shows  for  each  subaperture  the  phase  measurement 
(wavefront  slope  function) ,  displayed  in  the  form  of  3-D 
isometric  plots  and  contour  maps 

Figure  11  shows  the  effect  of  least  square  fitting  in  the 
overlap  areas.  The  horizontal  line  indicates  the  boundary 
between  two  successive  subapertures. 

Figure  12  shows  the  complete  synthesized  surface  (slope 
function)  displayed  in  the  form  of  3-D  isometric  plots  and 
contour  maps. 


6.4  Subaperture  Testing:  Final  Results 

The  FSIS  system  is  characterized  by  fully  automated  operation, 
whereby  successive  subaperture  inter ferograms  are  captured, 
computed  and  fitted  to  yield  a  full-surface,  multi-trace  wavefront 
measurement . 

The  raw  phase  measurement  yields  the  wavefront  slope  and  a 
numerical  integration  in  the  direction  of  the  long  side  of  the 
mirror  yield  a  set  of  actual  wavefront  profiles  or  traces. 
Interactive  processing  permits  to  extract  individual  traces  as 
desired. 

Figures  13  and  14  show  the  result  of  a  complete  measurement, 
yielding  about  100  traces  in  parallel. 

Note  that,  in  Figures  11  to  14,  the  continuity  of  the  contour  lines 
indicates  the  continuity  of  the  surface  and  thus  the  accuracy  of 
the  fit. 
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Figure  8:  Photographs  of  Interferograms  for  3  Subapertures 
FOR  Various  Alignment  Conditions 
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Subaperture  1:  Grating  in-Focus^  Mirror  Aligned 
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Subaperture  1:  Grating  in-Focus,  Mirror  Misaligned 
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Subaperture  1;  Grating  in-Focus^  Mirror  Strongly  Misaligned 


SlIBAPERTURE  ]:  GRATING  OUT-OF-FoCUS^  MIRROR  ALI 


SlIBAPERTLIRE  1:  GRATING  OUT-OF-FoCUS,  MIRROR  MISALIGNED  (POSITIVE) 
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Subaperture  2:  Grating  in-Focus,  Mirror  Aligned 
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SlIBArHRTURt  2:  GRATING  IN-F-QCIJS,  MiRROR  MISALIGNED 


Subaperture  2:  Grating  IN-Focusy  Mirror  Strongly  Misaligned 


Subaperture  3:  Grating  in-Focijs^  Mirror  Aligned 
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Figure  9:  Photographs  of  Subaperture  Interferograms  (5) 
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Subaperture  3 
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Subaperture  5 
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Figure  10;  Phase  Measurement  of  Individual  Subapertures 
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Phase  Measurement  of  Subaperture  Interferogram  #  1 
(3-D  Isometric  Plot) 


Phase  Measurement  of  Subaperture  Interferogram  #  1 
(Contour  Map) 


Phase  Measurement  of  Subaperture  Interferogram  #2 
(3-D  Isometric  Plot) 
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Phase  Measurement  of  Subaperture  Interferogram  §  2 
(Contour  Map) 


Phase  Measurement  of  Subaperture  Interferogram  #  3 
(3-D  Isometric  Plot) 


Phase  Measurement  of  Subaperture  Interferpgram  #3 
(Contour  Map) 


Phase  Measurement  of  Subaperture  Interferogram  U  4 
(3-D  Isometric  Plot) 
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Phase  Measurement  of  Subaperture  Interferogram  #  ^ 
(Contour  Map) 


Phase  Measurement  of  Subaperture  Interferogram  #  5 
(3-D  Isometric  Plot) 


-102- 


Phase  Measurement  of  Subaperture  Interferogram  #  5 
(Contour  Map) 
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Figure  11:  Phase  Measurement  Showing  Least  Square  Fitting 
IN  Subaperture  Overlap  Areas 
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Phase  Measurement  Showing  Least  Square  Fitting  in 
OVERLAP  Region  for  Subaperture  Interferograms  #  1  and  2 
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Phase  Measurement  Showing  Least  Square  Fitting  im 
Overlap  Region  for  Subaperture  Interferograms  #  2  and  3 
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Phase  Measurement  Showing  Least  Square  Fitting  in 
Overlap  Region  for  Subaperture  Interferograms  #  3  and  4 
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Figure  12:  Synthesis  of  Full  Surface  Phase  Measurement 
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Synthesis  of  Full-Surface  Phase  Measurement  from 
5  Subaperture  Interferograms  ^3-D  Isometric  Plot) 


Synthesis  of  Full-Surface  Phase  Measurement  from 
5  Subaperture  Interferograms  (Contour  Map) 
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Synthesis  of  Full-Surface  Phase  Measurement  from 
5  Subaperture  Interferograms  (Contour  Map) 
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Synthesis  of  Full-Surface  Phase  Measurement  from 
5  Subaperture  Interferograms  with  Tilt  Removed 
(Contour  Map) 
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Synthesis  of  Full-Surface  Phase  Measurement  from 
5  Subaperture  Interferograms  with  Tilt  Removed 
(Contour  Map) 


Figure  13:  Results  of  Subaperture  Testing  (Slope  Function) 
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LINDRICAL  MIRROR  :  MULTITRACE  MEASUREMENT  /  SLOPE  FUNCTION 
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CVLINDRICAL  MIRROR  :  MULT  I  TRACE  MEASUREMENT  /  SLOPE  FUNCTION 


-120 


HEftSUREHENT  /  SLOPE  FUNCTION 
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CVLINDRICAL  MIRROR  :  MULTITRACE  MEASUREMENT  /  SLOPE  FUNCTION 
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Figure  Results  of  Subaperture  Testing  (Wavefront  Function) 
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MULT  I TRACE  MEASUREMENT  /  WAVEFRONT  FUNCTION 


-141- 


142- 


NDRICAL  MIRROR  :  MULT  I  TRACE  MEASUREMENT  /  WAUEFRONT  FUN 
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LINDRICAL  MIRROR  :  MULT I TRACE  MEASUREMENT  /  WAVEFRONT  FUNCTION 
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CVLINDRICAL  MIRROR  :  MULT  I  TRACE  MEASUREMENT  /  WAUEFRONT  FUNCTION 
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CVLINDRICAL  MIRROR  :  MULT  I  TRACE  MEASUREMENT  /  UAUEFRONT  FUNCTION 
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